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cs| ! Abstract 

It is noted that one is now in possession of a set of facts, which may be viewed as 
■ the matching pieces of a puzzle ; in that all of them can be resolved by just one 

idea - that is grand unification. These include : (i) the observed family-structure, 
J> \ (ii) quantization of electric charge, (iii) meeting of the three gauge couplings, 

(iv) neutrino oscillations; in particular the mass of u-r (suggested by SuperK), 



(N 
00 



' (v) the intricate pattern of the masses and mixings of the fermions, including the 
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smallness of Vcb and the largeness of , and (vi) the need for B-L to implement 
baryogenesis (via leptogenesis) . All these pieces fit beautifully together within a 
\ single puzzle board framed by supersymmetric unification, based on SO(IO) or 

a string-unified G(224)-symmetry. The one and the most notable piece of the 
puzzle still missing, however, is proton decay. 

A concrete proposal is presented, within a predictive SO(10)/G(224)- 
. framework, that successfully describes the masses and mixings of all fermions, 

including the neutrinos - with eight predictions, all in agreement with observa- 
^ \ tion. Within this framework, a systematic study of proton decay is carried out, 

;h ' which pays special attention to its dependence on the fermion masses, including 

^ the superheavy Majorana masses of the right-handed neutrinos, and the thresh- 

old effects. The study (based on prior work and a recent update) shows that a 
conservative upper limit on the proton lifetime is about (1/2 - l)xlO^^ yrs, with 
vK^ being the dominant decay mode, and as a distinctive feature, fi^ being 
prominent. This in turn strongly suggests that an improvement in the current 
sensitivity by a factor of five to ten (compared to SuperK) ought to reveal proton 
decay. Otherwise some promising and remarkably successful ideas on unification 
would suffer a major setback. 



*Invited talk presented at the International School held at Erice Italy (Sept. 2000) and at the 
Dirac Medalists' symposium held at Trieste, Italy (Nov. 2000). 
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I. INTRODUCTION 



The standard model of particle physics, based on the gauge symmetry SU{2)l x U{1)y x 
SU{3)c is in excellent agreement with observations, at least up to energies of order 
100 GeV. Its success in turn constitutes a triumph of quantum field theory, especially of the 
notions of gauge invariance, spontaneous symmetry breaking, and renormalizability. The 
next step in the unification-ladder is associated with the concept of "grand unification", 
which proposes a unity of quarks and leptons, and simultaneously of their three basic forces: 
weak, electromagnetic and strong This concept was introduced on purely aesthetic 

grounds, in fact before any of the empirical successes of the standard model was in place. 
It was realized in 1972 that the standard model judged on aesthetic merits has some major 
shortcomings For example, it puts members of a family into five scattered multiplets, 

assigning rather peculiar hypercharge quantum numbers to each of them, without however 
providing a compelling reason for doing so. It also does not provide a fundamental reason 
for the quantization of electric charge, and it does not explain why the electron and proton 
possess exactly equal but opposite charges. Nor does it explain the co-existence of quarks 
and leptons, and that of the three gauge forces - weak, electromagnetic and strong - with 
their differing strengths. 

The idea of grand unification was postulated precisely to remove these shortcomings. It 
introduces the notion that quarks and leptons are members of one family, linked together by a 
symmetry group G, and that the weak, electromagnetic and strong interactions are aspects of 
one force, generated by gauging this symmetry G. The group G of course inevitably contains 
the standard model symmetry G(213) = SU{2)l x f/(l)y x SU{'i)c as a subgroup. Within 
this picture, the observed differences between quarks and leptons and those between the three 
gauge forces are assumed to be low-energy phenomena that arise through a spontaneous 
breaking of the unification symmetry G to the standard model symmetry G(213), at a very 
high energy scale M ^ ITeV. As a prediction of the hypothesis, such differences must then 
disappear and the true unity of quarks and leptons and of the three gauge forces should 
manifest at energies exceeding the scale M. 

The second and perhaps the most dramatic prediction of grand unification is proton 
decay. This important process, which would provide the window to view physics at truly 
short distances (< 10^'^° cm), is yet to be seen. Nevertheless, as I will stress in this talk, 
there has appeared over the years an impressive set of facts, favoring the hypothesis of grand 
unification. These include: 

(a) The observed family structure : The five scattered multiplets of the standard 
model, belonging to a family, neatly become parts of a whole {a single multiplet), with 
their weak hypercharges precisely predicted by grand unification. Realization of this feature 
calls for an extension of the standard model symmetry G(213) = SU(2)j;^xU(l)y xSU(3)'" 
minimally to the symmetry group G(224) = SU(2)iXSU(2)j:jxSU(4)'" which can be ex- 
tended further into the simple group SO(IO) 0|, but not SU(5) 0|. The G(224) symmetry 
in turn introduces some additional attractive features (see Sec.|l|), including especially the 
right-handed (RH) neutrinos (i'rs) accompanying the left-handed ones (i^l's), and B-L as a 
local symmetry. As we will see, both of these features now seem to be needed on empirical 
grounds. 
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(b) Meeting of the gauge couplings : Such a meeting is found to occur at a scale 
Mx ~ 2 X 10^^ GeV, when the three gauge couplings are extrapolated from their values 
measured at LEP to higher energies, in the context of supersymmetry |]^. This dramatic 
phenomenon supports the ideas of both grand unification and supersymmetry . These in 
turn may well emerge from a string theory |^ or M-theory (see discussion in Sec.[illD- 

(c) Mass of z/^ ~ 1/20 eV : Subject to the well-motivated assumption of hierar- 
chical neutrino masses, the recent discovery of atmospheric neutrino-oscillation at Su- 
per Kamiokande [0 suggests a value for m{vr) ~ l/2QeV. It has been argued (see e.g. 
Ref. that a mass of z/^ of this magnitude can be understood very simply by utilizing the 
SU(4)-color relation '^(i^T)i3i];ac ~ "^top ^^"^ SUSY unification scale Mx, noted above 
(See Sec.0). 

(d) Some intriguing features of fermion masses and mixings: These include: 
(i) the "observed" near equality of the masses of the b-quark and the r-lepton at the 
unification- scale (i.e. m° k, m°) and (ii) the observed largeness of the i^^-i^r oscillation 
angle (sin^ 2^°^^^^,^ > 0.83) together with the smallness of the corresponding quark mix- 
ing parameter V^f,(~ 0.04) ||TE[. As shown in recent work by Babu, Wilczek and me [|T3|, it 
turns out that these features and more can be understood remarkably well (see discussion 
in Sec.[V|) within an economical and predictive SO(10)-framework based on a minimal Higgs 
system. The success of this framework is in large part due simply to the group-structure of 
SO (10). For most purposes, that of G(224) suffices. 

(e) Baryogenesis : To implement baryogenesis successfully, in the presence of 
electroweak sphaleron effects [l^ , which wipe out any baryon excess generated at high tem- 
peratures in the (i?-L)-conserving mode, it has become apparent that one would need B-L as 
a generator of the underlying symmetry, whose spontaneous violation at high temperatures 
would yield, for example, lepton asymmetry (leptogenesis). The latter in turn is converted 
to baryon-excess at lower temperatures by electroweak sphalerons. This mechanism, it turns 
out, yields even quantitatively the right magnitude for baryon excess Jl^. The need for B-L, 
which is a generator of SU(4)-color, again points to the need for G(224) or SO(IO) as an 
effective symmetry near the unification-scale Mx- 

The success of each of these five features (a)-(e) seems to be non-trivial. Together they 
make a strong case for both supersymmetric grand unification and simultaneously for the 
G(224)/SO(10)-route to such unification, as being relevant to nature at short distances. 
However, despite these successes, as long as proton decay remains undiscovered, the hallmark 
of grand unification - that is quark-lepton trans formahility - would remain unrevealed. 

The relevant questions in this regard then are : What is the predicted range for the 
lifetime of the proton - in particular an upper limit - within the empirically favored route 
to unification mentioned above? What are the expected dominant decay modes within 
this route? Are these predictions compatible with current lower limits on proton lifetime 
mentioned above, and if so, can they still be tested at the existing or possible near-future 
detectors for proton decay? 

Fortunately, we are in a much better position to answer these questions now, compared 
to a few years ago, because meanwhile we have learnt more about the nature of grand 
unification. As noted above (see also Sec.^and Sec.|rV|), the neutrino masses and the meeting 
of the gauge couplings together seem to select out the supersymmetric G(224)/SO(10)-route 
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to higher unification. The main purpose of my talk here will therefore be to address the 
questions raised above, in the context of this route. For the sake of comparison, however, I 
will state the corresponding results for the case of supersymmetric SU(5) as well. 

My discussion will be based on a recent study of proton decay by Babu, Wilczek and me 
and an update of the same as presented here. Relative to other analysis, this study has 



three distinctive features: 

(a) It systematically takes into account the link that exists between proton decay and 
the masses and mixings of all fermions, including the neutrinos. 

(b) In particular, in addition to the contributions from the so-called "standard" d = 5 
operators [0 (see Sec.|VJ), it includes those from a new set of ci = 5 operators, related to 
the Majorana masses of the RH neutrinos |jl9[. These latter are found to be as important as 
the standard ones. 

(c) The work also incorporates GUT-scale threshold effects, which arise because of mass- 
splittings between the components of the SO(10)-multiplets, and lead to differences between 
the three gauge couplings. 

Each of these features turn out to be crucial to gaining a reliable insight into the nature 
of proton decay. Our study shows that the inverse decay rate for the VK~^-mode, which 
is dominant, is less than about 5 x 10^^ yrs for the case of MSSM embedded in SO(IO). 
This upper bound is obtained by making generous allowance for uncertainties in the matrix 
element and the SUSY-spectrum. Typically, the lifetime should of course be less than this 
bound. 

Proton decay is studied also for the case of the extended supersymmetric standard model 
(ESSM), that has been proposed a few years ago |^ on theoretical grounds, pertaining to the 
issues of string-unification and dilaton stabilization (see Sec.^ and the appendix). This case 
adds an extra pair of vector-like families at the TeV-scale, transforming as 16 + 16 of SO(IO), 
to the MSSM spectrum. While the case of ESSM is fully compatible with both neutrino- 
counting at LEP and precision electroweak tests, it can of course be tested directly at the 
LHC Our study shows that, with the inclusion of only the standard d=5 operators (defined 
in Sec.^lD, ESSM, embedded in SO(IO), can quite plausibly lead to proton lifetimes in the 
range of 10'^^ — 10^^ years, for nearly central values of the parameters pertaining to the SUSY- 
spectrum and the matrix element. Allowing for a wide variation of the parameters, owing 
to the contributions from both the standard and the neutrino mass-related d=5 operators 
(discussed in Sec.[VlD, proton lifetime still gets bounded above by about 10^^ years, even for 
the case of ESSM, embedded in SO (10) or a string - G(224). 

For either MSSM and ESSM, due to contributions from the new operators, the fi'^K^- 
mode is found to be prominent, with a branching ratio typically in the range of 10-50%. 
By contrast, minimal SUSY SU(5), for which the new operators are absent, would lead to 
branching ratios < 10~^ for this mode. 

Thus our study of proton decay, correlated with fermion masses, strongly suggests that 
discovery of proton decay should be imminent. In fact, one expects that at least candidate 
events should be observed in the near future already at SuperK. However, allowing for the 
possibility that the proton lifetime may well be closer to the upper bound stated above, a 
next-generation detector providing a net gain in sensitivity in proton decay-searches by a 
factor of 5-10, compared to SuperK, would certainly be needed not just to produce proton- 
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decay events, but also to clearly distinguish them from the background. It would of course 
also be essential to study the branching ratios of certain sub-dominant but crucial decay 
modes, such as the ^'^K^. The importance of such improved sensitivity, in the light of the 
successes of supersymmetric grand unification, is emphasized at the end. 



II. ADVANTAGES OF THE SYMMETRY G(224) AS A STEP TO HIGHER 

UNIFICATION 

As mentioned in the introduction, the hypothesis of grand unification was introduced to 
remove some of the conceptual shortcomings of the standard model (SM). To illustrate the 
advantages of an early suggestion in this regard, consider the five standard model multiplets 
belonging to the electron-family as shown : 

- — 1 

Here the superscripts denote the respective weak hypercharges IV (where Qem = -^3L+^iy/2) 
and the subscripts L and R denote the chiralities of the respective fields. If one asks : how 
one can put these five multiplets into just one multiplet, the answer turns out to be simple 
and unique. As mentioned in the introduction, the minimal extension of the SM symmetry 
G(213) needed, to achieve this goal, is given by the gauge symmetry : 

G(224) = SU(2)i X SU(2)k x SU(4)^. (2) 

Subject to left-right discrete symmetry {L ^ R), which is natural to G(224), all members 
of the electron family fall into the neat pattern : 



drp dii dh 



(3) 



The multiplets F| and F|. are left-right conjugates of each other and transform respectively 
as (2,1,4) and (1,2,4) of G(224); likewise for the muon and the tau families. Note that the 
symmetries SU(2)i and SU(2)ij are just like the familiar isospin symmetry, except that they 
operate on quarks and well as leptons, and distinguish between left and right chiralities. 
The left weak-isospin SU(2)i treats each column of F£ as a doublet; likewise SU(2)j^ for F^. 
The symmetry SU(4)-color treats each row of F£ and F|, as a quartet; thus lepton number is 
treated as the fourth color. Note also that postulating either SU(4)-color or SU(2)/j forces one 
to introduce a right-handed neutrino (uji) for each family as a singlet of the SM symmetry. 
This requires that there be sixteen two- component fermions in each family, as opposed to 
fifteen for the SM. The symmetry G(224) introduces an elegant charge formula : 



Qem = hi + ISR + ^ (4) 

expressed in terms of familiar quantum numbers I^l, hR and B-L, which applies to all forms 
of matter (including quarks and leptons of all six flavors, gauge and Higgs bosons). Note 
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that the weak hypercharge given by lV/2 = I^r + is now completely determined for 
all members of the family. The values of Yiy thus obtained precisely match the assignments 
shown in Eq. (p. Quite clearly, the charges I^l, I^ji and B-L, being generators respectively 
of 311(2)2,, SU(2)^ and SU(4)'^, are quantized; so also then is the electric charge Qem- 

In brief, the symmetry G(224) brings some attractive features to particle physics. These 
include : 

i) Unification of all 16 members of a family within one left-right self-conjugate multiplet; 

ii) Quantization of electric charge, with a reason for the fact that Qciectron = — Qproton 

iii) Quark- lepton unification (through SU(4) color); 

iv) Conservation of parity at a fundamental level P,pl[l; 

v) Right-handed neutrinos (i^^s) as a compelling feature; and 

vi) B-L local symmetry. 

As mentioned in the introduction, the two distinguishing features of G(224) - i.e. the exis- 
tence of the RH neutrinos and B-L as a local symmetry - now seem to be needed on empirical 
grounds. Furthermore, SU(4)-color provides simple relations between the masses of quarks 
and leptons, especially of those in the third family. As we will see in Secs.^ and 0, these 
are in good accord with observations. 

Believing in a complete unification, one is led to view the G(224) symmetry as part 
of a bigger symmetry, which itself may have its origin in an underlying theory, such as 
string theory. In this context, one may ask : Could the effective symmetry below the string 
scale in four dimensions (see Sec.fT|) be as small as just the SM symmetry G(213), even 
though the latter may have its origin in a bigger symmetry, which lives only in higher 
dimensions? I will argue in Sec.^ that the data on neutrino masses and the need for 
baryogenesis provide an answer to the contrary, suggesting that it is the effective symmetry 
in four dimensions, below the string scale, which must minimally contain either G(224) or a 
close relative G(214) = SU(2)Lxl3fixSU(4)^. 

One may also ask : does the effective four dimensional symmetry have to be any bigger 
than G(224) near the string scale? In preparation for an answer to this question, let us recall 
that the smallest simple group that contains the SM symmetry G(213) is SU(5) [Q. It has 
the virtue of demonstrating how the main ideas of grand unification, including unification 
of the gauge couplings, can be realized. However, SU(5) does not contain G(224) as a 
subgroup. As such, it does not possess some of the advantages listed above. In particular, 
it does not contain the RH neutrinos as a compelling feature, and B-L as a local symmetry. 
Furthermore, it splits members of a family into two multiplets : 5-1-10. 

By contrast, the symmetry SO(IO) has the merit, relative to SU(5), that it contains 
G(224) as a subgroup, and thereby retains all the advantages of G(224) listed above. (As a 
historical note, it is worth mentioning that these advantages had been motivated on aesthetic 
grounds through the symmetry G(224) and all the ideas of higher unification were in 
place before it was noted that G (224) (isomorphic to SO(4)xSO(6)) embeds nicely into 
SO(IO) 0). Now, SO(IO) even preserves the 16-plet family- structure of G(224) without a 
need for any extension. By contrast, if one extends G(224) to the still higher symmetry 
Eg |2^, the advantages (i)-(vi) are retained, but in this case, one must extend the family- 
structure from a 16 to a 27-plet, by postulating additional fermions. In this sense, there 
seems to be some advantage in having the effective symmetry below the string scale to be 
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minimally G(224) (or G(214)) and maximally no more than SO(IO). I will compare the 
relative advantage of having either a string-derived G(224) or a string-SO(lO), in the next 
section. First, I discuss the implications of the data on coupling unification. 



III. THE NEED FOR SUPERSYMMETRY : MSSM VERSUS STRING 

UNIFICATIONS 

It has been known for some time that the precision measurements of the standard model 
coupling constants (in particular sin^ 9w) at LEP put severe constraints on the idea of 
grand unification. Owing to these constraints, the non-supersymmetric minimal SU(5), and 
for similar reasons, the one-step breaking minimal non-supersymmetric SO(10)-model as 



well, are now excluded But the situation changes radically if one assumes that the 

standard model is replaced by the minimal supersymmetric standard model (MSSM), above 
a threshold of about 1 TeV. In this case, the three gauge couplings are found to meet 0], to 
a very good approximation, barring a few percent discrepancy which can be attributed to 
threshold corrections (see Appendix). Their scale of meeting is given by 

Mx^2x 10^^ GeV (MSSM or SUSY SU(5)) (5) 

This dramatic meeting of the three gauge couplings, or equivalently the agreement of the 
MSSM-based prediction of sin^ Ow{^z)'y\i = 0.2315 ± 0.003 |2^ with the observed value of 



sin Owimz) = 0.23124 ± 0.00017 [T^, provides a strong support for the ideas of both grand 
unification and supersymmetry, as being relevant to physics at short distances. 

In addition to being needed for achieving coupling unification there is of course an inde- 
pendent motivation for low-energy supersymmetry - i.e. for the existence of SUSY partners 
of the standard model particles with masses of order 1 TeV. This is because it protects the 
Higgs boson mass from getting large quantum corrections, which would (otherwise) arise 
from grand unification and Planck scale physics. It thereby provides at least a technical 
resolution of the so-called gauge-hierarchy problem. In this sense low-energy supersymmetry 
seems to be needed for the consistency of the hypothesis of grand unification. Supersym- 
metry is of course also needed for the consistency of string theory. And most important, 
low-energy supersymmetry can be tested at the LHC, and possibly at the Tevatron. 

The most straightforward interpretation of the observed meeting of the three gauge cou- 
plings and of the scale Mx, is that a supersymmetric grand unification symmetry (often 
called GUT symmetry), like SU(5) or SO(IO), breaks spontaneously at Mx into the stan- 
dard model symmetry G(213). 

Even if supersymmetric grand unification may well be a good effective theory below a 
certain scale M > Mx, it ought to have its origin within an underlying theory like string/M 
theory. Such a theory is needed to unify all the forces of nature including gravity, and to 
provide a good quantum theory of gravity. It is also needed to provide a rationale for the 
existence of flavor symmetries (not available within grand unification), which distinguish 
between the three families and can resolve certain naturalness problems including those 
associated with inter-family mass hierarchy. 

In the context of string or M theory, an alternative interpretation of the observed meeting 
of the gauge couplings is however possible. This is because, even if the effective symmetry in 
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four dimensions emerging from a higher dimensional string theory is non-simple, like G(224) 
or G(213), string theory can still ensure familiar unification of the gauge couplings at the 
string scale. In this case, however, one needs to account for the small mismatch between the 
MSSM unification scale Mx (given above), and the string unification scale, given by Mst ~ 
Qst X 5.2 X 10^^ GeV ^ 3.6 x 10^^ GeV (Here we have put a^t = aGuAMSSM) ^ 0.04) [^. 
Possible resolutions of this mismatch have been proposed. These include : (i) utilizing the 
idea of string- duality |^6[ which allows a lowering of Mg^ compared to the value shown above, 
or alternatively (ii) the idea of a semi-perturbative unification that assumes the existence of 
two vector-like families, transforming as (16 + 16) of SO(IO), with masses of order one TeV 
pO| . The latter raises acuT to about 0.25-0.3 and simultaneously Mx, in two loop, to about 
(1/2 — 2) X 10^^ GeV. (Other mechanisms resolving the mismatch are reviewed in Ref. ||27|| ). 
In practice, a combination of the two mechanisms mentioned above may well be relevant. Q 

While the mismatch can thus quite plausibly be removed for a non-GUT string-derived 
symmetry like G(224) or G(213), a GUT symmetry like SU(5) or SO(IO) would have an 
advantage in this regard because it would keep the gauge couplings together between Mgt and 
Mx (even if Mx ~ Mgt /20), and thus not even encounter the problem of a mismatch between 
the two scales. A supersymmetric GUT-solution (like SU(5) or SO(IO)), however, has a 
possible disadvantage as well, because it needs certain color triplets to become superheavy 
by the so-called doublet-triplet sphtting mechanism (see Sec.|V5 and Appendix), in order to 
avoid the problem of rapid proton decay. However, no such mechanism has emerged yet, in 
string theory, for the GUT-like solutions PB| . 

Non-GUT string solutions, based on symmetries like G(224) or G(2113) for example, have 
a distinct advantage in this regard, in that the dangerous color triplets, which would induce 
rapid proton decay, are often naturally projected out for such solutions Furthermore, 



the non-GUT solutions invariably possess new "flavor" gauge symmetries, which distinguish 
between families. These symmetries are immensely helpful in explaining qualitatively the 
observed fermion mass-hierarchy (see e.g. Ref. [Q) and resolving the so-called naturalness 



problems of supersymmetry such as those pertaining to the issues of squark-degeneracy [31 



CP violation |^ and quantum gravity-induced rapid proton decay ||3c 

Weighing the advantages and possible disadvantages of both, it seems hard at present 
to make a priori a clear choice between a GUT versus a non-GUT string-solution. As 
expressed elsewhere [Q, it therefore seems prudent to keep both options open and pursue 
their phenomenological consequences. Given the advantages of G(224) or SO(IO) in the light 



^ I have in mind the possibility of string-duality ||2^ lowering Mst for the case of semi-perturbative 
unification (for which agt ~0.25, and thus, without the use of string-duality, Mst would be about 
10^^ GeV) to a value of about (l-2)xl0^^ GeV (say), and semi-perturbative unification |]2^ raising 
the MSSM value of Mx to about 5x10^*^ GeV« M^t(l/2 to 1/4) (say). In this case, an intermediate 
symmetry like G(224) emerging at Mst would be effective only within the short gap between Mst 
and Mx, where it would break into G(213). Despite this short gap, one would still have the benefits 
of SU(4)-color that are needed to understand neutrino masses (see sec. 4). At the same time, since 
the gap is so small, the couplings of G(224), unified at Mgt would remain essentially so at Mx, so 
as to match with the "observed" coupling unification, of the type suggested in Ref. |2C 
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of the neutrino masses (see Secs.^ and [TVD , I will thus proceed by assuming that either a 
suitable G(224)-solution with a mechanism of the sort mentioned above, or a realistic SO(IO)- 
solution with the needed doublet-triplet mechanism, will emerge from string theory. We 
will see that with this broad assumption, an economical and predictive framework emerges, 
which successfully accounts for a host of observed phenomena, and makes some crucial 
testable predictions. Fortunately, it will turn out that there are many similarities between 
the predictions of a string-unified G(224) and SO(IO) frameworks, not only for the neutrino 
and the charged fermion masses, but also for proton decay. I next discuss the implications 
of the mass of suggested by the SuperK data. 



IV. MASS OF Vr-. EVIDENCE IN FAVOR OF THE G(224) ROUTE 



One can obtain an estimate for the mass of v}^ in the context of G(224) or SO(IO) by 
using the following three steps (see e.g.Ref. |T^): 



(i) Assume that B— L and /s/j, contained in a string-derived 0(224) or SO(IO), break 
near the unification-scale: 

Mx ~ 2 X 10^^ OeV , (6) 

through VEVs of Higgs multiplets of the type suggested by string-solutions - i.e. ((1, 2, A)h) 
for 0(224) or (16/^) for SO(IO), as opposed to 126// which seems to be unobtainable (at 
least) in weakly interacting string theory |3^. In the process, the RH neutrinos (z^lj), which 



are singlets of the standard model, can and generically will acquire superheavy Majorana 
masses of the type M]j v^ji uj^, by utilizing the VEV of (16//) and effective couplings of 
the form: 



Cm {SO{10)) = fij 16, ■ 16, 16// • 16///M + h.c. (7) 

A similar expression holds for 0(224). Here i,j = 1,2,3, correspond respectively to 
e, fi and r families. Such gauge-invariant non-renormalizable couplings might be expected 
to be induced by Planck-scale physics, involving quantum gravity or stringy effects and/or 
tree-level exchange of superheavy states, such as those in the string tower. With fij (at 
least the largest among them) being of order unity, we would thus expect M to lie between 
Mpianck ~ 2 X 10^*^ OeV and Mstring ~ 4 X 10^'' OeV. Ignoring for the present off-diagonal 
mixings (for simplicity), one thus obtains 0: 

M^R ^ ^^^J^ ^ /33 (2 X 10^^ OeV) p2 (Mpianck/M) (8) 

This is the Majorana mass of the RH tau neturino. Ouided by the value of Mx, we have 
substituted (l6//) = (2 x 10^^ OeV) p ,with p ^ 1/2 to 2(say). 



^ The effects of neutrino-mixing and of possible choice of M = Mstring ~ 4 x 10^^ GeV (instead 
of M = M Planck) on M^p are considered in Ref. |Q. 
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(ii) Now using SU(4)-color and the Higgs multiplet (2, 2, 1)h of G(224) or equivalently 
IOh of SO(IO), one obtains the relation mr{Mx) = mh{Mx), which is known to be successful. 
Thus, there is a good reason to believe that the third family gets its masses primarily from 
the IOh or equivalently (2,2, 1)^^ (see sec. 5). In turn, this implies: 



Dirac) 



mtop(Mx) ^ (100-120)GeV (9) 



Note that this relationship between the Dirac mass of the tau-neutrino and the top-mass is 
special to SU(4)-color. It does not emerge in SU(5). 

(iii) Given the superheavy Majorana masses of the RH neutrinos as well as the Dirac 
masses as above, the see-saw mechanism yields naturally light masses for the LH neu- 
trinos. For 1^1 (ignoring flavor- mixing) , one thus obtains, using Eqs.(^ and (^, 



\2 



^ ^ [(l/20)eV(l-1.44)//33p'](M/M„,„,,) (10) 

Now, assuming the hierarchical pattern ^«(z/£) <^ ^(^i^) -C m{i''£), which is suggested by 
the see-saw mechanism, and further that the SuperK observation represents z/^ — uj^ (rather 
than z/{^ — ux) oscillation, the observed ^ 1/2(10"^- 10~^) eV^ corresponds to m{i'1) ~ 
(1/15 - 1/40) eV. It seems truly remarkable that the expected magnitude of m(z/2), given by 
Eq.(PDD, is just about what is suggested by the SuperK data, if fzj, P^i.^Pianck/ ~ 1-3 to 
1/2. Such a range for p^{Mpianck/M) seems most plausible and natural (see discussion in 
Ref. [0). Note that the estimate ([l0| ) crucially depends upon the supersymmetric unification 
scale, which provides a value for M^r, as well as on SU(4)-color that yields m{i''^^^^^.The 
agreement between the expected and the SuperK results thus clearly favors supersymmetric 
unification, and in the string theory context, it suggests that the effective symmetry below the 
string-scale should contain SU(4)-color. Thus, minimally this effective symmetry should be 
either G(214) or G(224), and maximally as big as SO(IO), if not Eg. 

By contrast, if SU(5) is regarded as either a fundamental symmetry or as the effective 
symmetry below the string scale, there would be no compelling reason based on symmetry 
alone, to introduce a ur, because it is a singlet of SU(5). Second, even if one did introduce 
z/]j by hand, their Dirac masses, arising from the coupling 5i{5H)i^R, would be unrelated 
to the up- flavor masses and thus rather arbitrary (contrast with Eq. (^). So also would be 
the Majorana masses of the i'r's, which are SU(5)-invariant, and thus can be even of order 
string scale . This would give m{v'^) in gross conflict with the observed value. 

Before passing to the next section, it is worth noting that the mass of v^- suggested by 
SuperK, as well as the observed value of sin^ (see SecPTTD, provide valuable insight into 
the nature of GUT symmetry breaking. They both favor the case of a single-step breaking 
(SSB) of SO(IO) or a string-unified G(224) symmetry at a scale of order Mx, into the 
standard model symmetry G(213), as opposed to that of a multi-step breaking (MSB). The 
latter would correspond, for example, to SO(IO) (or G(224)) breaking at a scale Mi into 
G(2213), which in turn breaks at a scale M2 « Mi into G(213). One reason why the 
case of single-step breaking is favored over that of multi-step breaking is that the latter 
can accommodate but not really predict sin^ 6'vy, whereas the former predicts the same 
successfully. Furthermore, since the Majorana mass of z/]j arises arises only after B — L and 
I^R break, it would be given, for the case of MSB, by M^r ~ /33(M|/M), where M ~ Mst 
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(say). If M2 < Mx ~ 2 X 10^^ GeV, and M> Mx, one would obtain too low a value (<< 10^^ 
GeV) for M^n (compare with Eq.(8)), and thereby too large a value for m(z/£), compared 
to that suggested by SuperK. By contrast, the case of SSB yields the right magnitude for 
mifr) (see Eq. (10)). 

Thus the success of the result on m(z/^) discussed above not only favors the symmetry 
G(224) or SO(IO), but also clearly suggests that B — L and I^r break near the conventional 
GUT scale Mx ~ 2 x 10^^ GeV, rather than at an intermediate scale << Mx- In other words, 
the observed values of both sin^ 6w and m^Ur) favor only the simplest pattern of symmetry- 
breaking, for which SO(IO) or a string-derived G(224) symmetry breaks in one step to the 
standard model symmetry, rather than in multiple steps. It is of course only this simple 
pattern of symmetry breaking that would be rather restrictive as regards its predictions for 
proton decay (to be discussed in Sec.[VT|). I next discuss the problem of understanding the 
masses and mixings of all fermions. 



V. UNDERSTANDING FERMION MASSES AND NEUTRINO OSCILLATIONS 

IN SO(IO) 

Understanding the masses and mixings of all quarks and charged leptons, in conjunction 
with those of the neutrinos, is a goal worth achieving by itself. It also turns out to be 
essential for the study of proton decay. I therefore present first a recent attempt in this 
direction, which seems most promising |T^. A few guidelines would prove to be helpful in 



this regard. The first of these is motivated by the desire for economy and the rest by data. 

1) Hierarchy Through Off-diagonal Mixings : Recall earlier attempts that 
attribute hierarchical masses of the first two families to mass matrices of the form : 

M=(° ;)mf), (11) 

for the {d,s) quarks, and likewise for the {u,c) quarks. Here e ~ 1/10. The hierarchical 
patterns in Eq. (0) can be ensured by imposing a suitable flavor symmetry which distin- 
guishes between the two families (that in turn may have its origin in string theory (see e.g. 



Ref I^O])- Such a pattern has the virtues that (a) it yields a hierarchy that is much larger 
than the input parameter e : (m^/ms) ~ <^ e, and (b) it leads to an expression for the 
cabibbo angle : 

e. 



c ^ 



(12) 



which is rather successful. Using ^Jm^jms 0.22 and \fm^Jrnc ~ 0.06, we see that Eq. 
(P^) works to within about 25% for any value of the phase 0. Note that the square root 
formula (like \fm^m'^ for the relevant mixing angle arises because of the symmetric form 
of M in Eq. (pT]), which in turn is ensured if the contributing Higgs is a 10 of SO(IO). A 
generalization of the pattern in Eq. ([TT| ) would suggest that the first two families (i.e. the 
e and the \x) receive masses primarily through their mixing with the third family (r), with 
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(1, 3) and (1, 2) elements being smaller than the (2, 3); while (2, 3) is smaller than the (3, 3). 
We will follow this guideline, except for the modification noted below. 

2) The Need for an Antisymmetric Component : Although the symmetric hier- 
archical matrix in Eq. (|Tl|) works well for the first two families, a matrix of the same form 
fails altogether to reproduce Vcb, for which it yields : 




m„ • irrir 
- 

nib V 




(13) 



Given that \fni^[mb ~ 0.17 and ^Jlnjnit ~ 0.0.06, we see that Eq. (|T3D would yield Vcb 
varying between 0.11 and 0.23, depending upon the phase x- This is too big, compared to 
the observed value of Vcb ~ 0.04 ±0.003, by at least a factor of 3. We interpret this failure as 
a clue to the presence of an antisymmetric component in M, together with symmetrical ones 

(so that ^ rriji), which would modify the relevant mixing angle to ^/——^/ — where 



rrij y mj. 



rrii and nij denote the respective eigenvalues. 

3) The Need for a Contribution Proportional to B-L : The success of the relations 
m° ~ m°, and m° ~ ''^{^T)Dirac (^^e Sec.|VD, suggests that the members of the third family 
get their masses primarily from the VEV of a SU(4)-color singlet Higgs field that is indepen- 
dent of B-L. This is in fact ensured if the Higgs is a 10 of SO(IO). However, the empirical 
observations of m° ~ m^/3 and m° ~ 3m^ clearly call for a contribution proportional to 



B-L as well. Further, one can in fact argue that the suppression of Vcb (in the quark-sector) 
together with an enhancement of 0'^'^^^,^^ (in the lepton sector) calls for a contribution that 
is not only proportional to B-L, but also antisymmetric in the family space (as suggested 
above in item (^). We show below how both of these requirements can be met, rather easily, 
in SO(IO), even for a minimal Higgs system. 

4) Up-Down Asymmetry: Finally, the up and the down-sector mass matrices must 
not be proportional to each other, as otherwise the CKM angles would all vanish. Note that 
the cubic couplings of a single IOh will not serve the purpose in this regard. 

Following Ref. []T1|, I now present a simple and predictive mass-matrix, based on SO(IO), 



that satisfies all four requirements (0), (0), (|) and (|). The interesting point is that one 
can obtain such a mass-matrix for the fermions by utilizing only the minimal Higgs system, 
that is needed anyway to break the gauge symmetry SO(IO). It consists of the set : 

Hminimal = {45//, 16//, 16//, 10//} . (14) 

Of these, the VEV of (45//) ~ Mx breaks SO(IO) into G(2213), and those of (16//) = 
(16//) ~ Mx break G(2213) to G(213), at the unification-scale Mx- Now G(213) breaks at 
the electroweak scale by the VEV of (10//) to U(l)e^x SU(3)^ 

One might have introduced large-dimensional tensorial multiplets of SO (10) like 126// and 
120//, both of which possess cubic level Yukawa couplings with the fermions. In particular, 
the coupling 16il6j(120//) would give the desired family-antisymmetric as well as [B-L)- 
dependent contribution. We do not however introduce these multiplets in part because 



they do not seem to arise in string solutions [35|, and in part also because mass-splittings 



within such large- dimensional multiplets could give excessive threshold corrections to 03 (m^ 
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(typically exceeding 20%), rendering observed coupling unification fortuitous. By contrast, 
the multiplets in the minimal set (shown above) do arise in string solutions leading to 
SO(IO). Furthermore, the threshold corrections for the minimal set are found to be naturally 
small, and even to have the right sign, to go with the observed coupling unification [|14| (see 
Appendix). 

The question is: can the minimal set of Higgs multiplets (see Eq.(|T4D) meet all the 
requirements listed above? Now IOh (even several lO's) can not meet the requirements of 
antisymmetry and (i?-L)-dependence. Furthermore, a single IOh cannot generate CKM- 
mixings. This impasse disappears, however, as soon as one allows for not only cubic, but 
also effective non-renormalizable quartic couplings of the minimal set of Higgs fields with the 
fermions. These latter couplings could of course well arise through exchanges of superheavy 
states (e.g. those in the string tower) involving renormalizable couplings, and/or through 
quantum gravity. 

Allowing for such cubic and quartic couplings and adopting the guideline (^ of hierar- 
chical Yukawa couplings, as well as that of economy, we are led to suggest the following 
effective lagrangian for generating Dirac masses and mixings of the three families |]T^ (for a 
related but different pattern, involving a non-minimal Higgs system, see Ref ||39|| ). 



^Yuk = h3 I63 I63 IOh + [ I62 I63 IOh + 033 I62 I63 IOh 45h/M 
+ c/23 I62 I63 16h 16h/M] + {ai2 I61 162 IOh 45h/M 

+ (712 I61 162 16h 16h/M} . (15) 

Here, M could plausibly be of order string scale. Note that a mass matrix having essentially 
the form of Eq. (pH]) results if the first term /i33(10h) is dominant. This ensures m° ^ 
m° and m° ~ m{v£,irac)^- Following the assumption of progressive hierarchy (equivalently 
appropriate flavor symmetries Q ),we presume that /123 ~ /133/IO, while /i22 and hu, which are 
not shown, are assumed to be progressively much smaller than /i23. Since (45//) ~ (16//) ~ 
Mx-i while M ~ Mst ~ lOMx, the terms 023 (45//) /M and (723 (16//) /M can quite plausibly 
be of order /133/IO, if 023 ~ (723 ~ ^33- By the assumption of hierarchy, we presume that 
ai2 < ^23, and gi2 < 5-23 

It is interesting to observe the symmetry properties of the 023 and (723-terms. Although 
10// X 45// = 10 + 120 + 320, given that (45//) is along B-L, which is needed to implement 



^Although no explicit string solution with the hierarchy in all the Yukawa couplings in Eq.([l5|) - 
i.e. in /ijj, aij and gij - exists as yet, one can postulate flavor symmetries of the type alluded to 
(e.g. two abelian U(l) symmetries), which assign flavor charges not only to the fermion families 
and the Higgs multiplets, but also to a few (postulated) SM singlets that acquire VEVs of order 
Mx. The flavor symmetry - allowed effective couplings such as I62I63IO// < S > /M would lead to 
^23 ~< S > /M ~ 1/10. One can verify that the full set of hierarchical couplings shown in Eq.(p!5[) 
can in fact arise in the presence of two such U(l) symmetries. String theory (at least) offers the 
scope (as indicated by the solutions of Refs. [^] and |2^) for providing a rationale for the existence 
of such flavor symmetries, together with that of the SM singlets. For example, there exist solutions 
with the top Yukawa coupling being leading and others being hierarchical (as in Ref. pO|). 
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doublet-triplet splitting (see Appendix), only 120 in the decomposition contributes to the 
mass-matrices. This contribution is, however, antisymmetric in the family-index and, at the 
same time, proportional to B-L. Thus the 023 term fulfills the requirements of both anti- 
symmetry and (B-L)- dependence, simultaneously^ . With only hij and ajj-terms, however, 
the up and down quark mass-matrices will be proportional to each other, which would yield 
VcKM = 1- This is remedied by the gij coupling, because, the 16h can have a VEV not only 
along its SM singlet component (transforming as z/^j) which is of GUT-scale, but also along 
its electroweak doublet component - call it 16^ - of the electroweak scale. The latter can 
arise by the the mixing of 16^ with the corresponding doublet (call it 10^) in the IOh- The 
MSSM doublet H^, which is light, is then a mixture of 10^ and 16^, while the orthogonal 
combination is superheavy (see Appendix). Since (16^) contributes only to the down- flavor 
mass matrices, but not to the up-flavor, the g23 and gu couplings generate non-trivial CKM- 
mixings. We thus see that the minimal Higgs system (as shown in Eg. (fl^j ) satisfies apriori 
all the qualitative requirements (l)-(4), including the condition ofVcKM 7^ 1- I now discuss 
that this system works well even quantitatively. 

With these six effective Yukawa couplings, the Dirac mass matrices of quarks and leptons 
of the three families at the unification scale take the form : 

U = \ -e' e + o- mjj, D = \ - 1' + 7]' e + r] m^, 
\0-e + al/ \0 -e + ri 1 J 



/ - 3e' \ / -3e' + T]' \ 

N = i 3e' -3e + (T mc/, L = 3e' + V -3e + \ tud. (16) 
\0 3e + (T 1 / \0 3e + ri 1 / 

Here the matrices are multiplied by left-handed fermion fields from the left and by anti- 
fermion fields from the right. {U, D) stand for the mass matrices of up and down quarks, 
while (A^, L) are the Dirac mass matrices of the neutrinos and the charged leptons. The 
entries l,e,and a arise respectively from the /i33,a23 and /123 terms in Eq. (|15|) , while rj 
entering into D and L receives contributions from both (723 and h23; thus i] ^ a. Similarly rj' 
and e' arise from gi2 and 012 terms respectively. Note the quark-lepton correlations between 
U and A^ as well as D and L, and the up-down correlations between U and D as well as A^ 
and L. These correlations arise because of the symmetry property of G(224). The relative 
factor of —3 between quarks and leptons involving the e entry reflects the fact that (45h) ~ 
to(B-L), while the antisymmetry in this entry arises from the group structure of SO(IO), as 
explained above^. As we will see, this e-entry helps to account for (a) the differences between 
nis and m^, (b) that between nid and rrie, and also, (c) the suppression of V^b together with 
the enhancement of the z/^-z/,- oscillation angle. 



^The analog of 10^-45// for the case of G(224) would be xh = (2, 2, 1)^ • (1, 1, 15)h- Although in 
general, the coupling of xh to the fermions need not be antisymmetric, for a string-derived G(224), 
the multiplet (1,1, 15)/^ is most Hkely to arise from an underlying 45 of SO(IO) (rather than 210); 
in this case, the couplings of xh must be antisymmetric like that of IOh ■ 45i^. 
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The mass matrices in Eq.(|T6|) contain 7 parameters 0: e, a, r], uid = h^^, (10^), mu = 
^33(10[/), y' and e'. These may be determined by using, for example, the following input 

1.37 GeV, m,(l GeV) = 110-116 MeV [|g, m„(l GeV) 



values: mS'^^'^ = 174 GeV, 



mAnir 



6 MeV and the observed masses of e, jj, and r, which lead to (see Ref. [Q, for details): 



a ~ 0.110 , r/ ~ 0.151 , e ~ - 0.095 , \i\ ^ 4.4 x 10"^ and e' ^ 2 x 10" 



m(7 



(100-120) GeV, mD ^ ^^(Mf/) ~ 1.5 GeV. 



;i7) 



Here, I will assume, only for the sake of simplicity, as in Ref. |Tj|, that the parameters 
are real 0. Note that in accord with our general expectations discussed above, each of the 
parameters cr, t] and e are found to be of order 1/10, as opposed to being [] 0(1) or O(10~^), 
compared to the leading (3,3)-element in Eq. ([I6|). Having determined these parameters, 
we are led to a total of five predictions involving only the quarks (those for the leptons are 
listed separately) : 



mi. 



m°(l - Se^); thus mb(mb) ~ (4.6-4.9) GeV 



IK 



cb\ 



\a — ■r]\ 



V 



1/2 



T] — e 
md(lGeV) ~ 8 MeV 



a 



a 



1/2 



0.045 



iKb/Kfcl - Vmjm'c ^ 0.07. 



(19) 
(20) 
(21) 
(22) 



In making these predictions, we have extrapolated the GUT-scale values down to low energies 
using a3{mz) = 0.118, a SUSY threshold of 500 GeV and tan/5 = 5. The results depend 
weakly on these choices, assuming tan/3 ~ 2-30. Further, the Dirac masses and mixings of 
the neutrinos and the mixings of the charged leptons also get determined. We obtain : 



^ Of these, ~ can in fact be estimated to within 20% accuracy by either using the argument 
of radiative electroweak symmetry breaking, or some promising string solutions (see e.g. Ref. [^0| ). 



^ Babu and I have recently studied supersymmetric CP violation within the G(224)/SO(10) frame- 
work, by using precisely the fermion mass-matrices as in Eq.(ll). We have observed (H that com- 
plexification of the parameters can lead to observed CP violation, without upsetting in the least 
the success of Ref. (i.e. of the fermion mass-matrices of Eq.(^)) in describing the masses and 
mixings of all fermions, including neutrinos. Even with complexification the relative signs and the 
approximate magnitudes of the real parts of the parameters must be the same as in Eq.([T7|), to 
retain the success. 



^ This is one characteristic difference between our work and that of Ref. |39|, where the (2,3)- 
element is even bigger than the (3,3). 
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m^^{Mu) ^ 100-120 GeV; m^^{Mu) - 8 GeV, (23) 

1/2 

~ 0.437 (24) 
m^^ ~ [9e'V(9e2 - (T^)]mu ~ 0.4 MeV (25) 

\J~mJrn^ ~ 0.85 \J~rnJrn^j_ ~ 0.06 (26) 

C - ^ \/"^e/m, (m^/m,) ~ 0.0012 . (27) 

In evaluating 6*^^, we have assumed e' and 77' to be relatively positive. 

Given the bizarre pattern of quark and lepton masses and mixings, it seems remark- 
able that the simple pattern of fermion mass-matrices, motivated by the group theory of 
G(224)/SO(10), gives an overall fit to all of them (Eqs.([T8|) through (^)) which is good 
to within 10%. This includes the two successful predictions on rrih and Vcb (Eqs.(|18|) and 
(p!9|)). Note that in supersymmetric unified theories, the "observed" value of mb{mi)) and 
renormalizat ion-group studies suggest that, for a wide range of the parameter tan/5, m° 
should in fact be about 10-20% lower than |^ . This is neatly explained by the rela- 
tion: m° ~ m°(l — 8e^) (Eq. (|TH|)), where exact equality holds in the limit e — > (due to 
SU(4)-color), while the decrease of m° compared to m° by 8e^ ~ 10% is precisely because 
the off-diagonal e-entry is proportional to B-L (see Eq. (|16D). 

Specially intriguing is the result on Vcb ~ 0.045 which compares well with the observed 
value of ~ 0.04. The suppression of Vcb, compared to the value of 0.17 ± 0.06 obtained from 
Eq. (|TB[), is now possible because the mass matrices (Eq. (|TBp) contain an antisymmetric 
component oc e. That corrects the square-root formula Ogb = \/mJmb (appropriate for 
symmetric matrices, see Eq. (pH])) by the asymmetry factor Kr] + / {1] — e)]^^"^ (see Eq. 
(19)), and similarly for the angle 9^- This factor suppresses Vcb if V ^ have opposite 
signs. The interesting point is that, the same feature necessarily enhances the corresponding 
mixing angle 9^^^ in the leptonic sector, since the asymmetry factor in this case is given by 
[(— 3e + r])/(3e + r/)]^''^ (see Eq. (24)). This enhancement of 6*^^ helps to account for the 
nearly maximal oscillation angle observed at SuperK (as discussed below). This intriguing 
correlation between the mixing angles in the quark versus leptonic sectors - that is suppres- 
sion of one implying enhancement of the other - has become possible only because of the 
e-contribution, which is simultaneously antisymmetric and is proportional to B-L. That in 
turn becomes possible because of the group-property of SO(IO) or a string-derived G(224)^. 

Taking stock, we see an overwhelming set of facts in favor of B-L and in fact for the 
full SU(4)-color-symmetry. These include: (i) the suppression of Vcb, together with the 
enhancement of just mentioned above, (ii) the successful relation ^ m°(l — 8e^), (iii) 
the usefulness again of the SU(4)-color-relation 'm{i>'^^^^J^ ~ m° in accounting for m(z/J^)( 
see Sec. 4 ), and (iv) the agreement of the relation |m°/m°| = |(e^ — ?7^)/(9e^ — ?7^)| with 
the data, in that the ratio is naturally less than 1, ii rj ^ e. The presence of 9e^ in the 
denominator is because the off-diagonal entry is proportional to B-L. Finally, the need for 
(B-L)- as a local symmetry, to implement baryogenesis, has been noted in Sec.l. 



9' ^ 



3e + 1] ^ \J m^/rrij 



- 3e + ?7 



9' ~ 



7] 



3e' 



T]' + 3e' 
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Turning to neutrino masses, while all the entries in the Dirac mass matrix are now 
fixed, to obtain the parameters for the light neutrinos, one needs to specify those of the 
Majorana mass matrix of the RH neutrinos (z/^''^'^). Guided by economy and the assumption 
of hierarchy, we consider the following pattern : 




M^ = [ y \ Mr. (28) 
As discussed in Sec.^, the magnitude of Mr ~ (5-15) x 10^^ GeV can quite plausibly 



be justified in the context of supersymmetric unificatonQ (e.g. by using M ^ Mgt ~ 4 x 10^"^ 
GeV in Eq. (||)). To the same extent, the magnitude of miyr) ~ (1/10-1/30) eV, which is 
consistent with the SuperK value, can also be anticipated. Thus there are effectively three 
new parameters: x, and z. Since there are six observables for the three light neutrinos, 
one can expect three predictions. These may be taken to be 6*°^^^, m^^ (see Eq. (|10D), and 
for example 

Assuming successively hierarchical entries as for the Dirac mass matrices, we presume 
that \y\ ~ 1/10, \z\ < |i/|/10 and |x| < z"^ . Now given that m{yr) ~ 1/20 eV (as estimated in 
Eq. (|ToD ), the MSW solution for the solar neutrino puzzle suggests that m(i/^)/m(z/T-) ^ 
1/10-1/30. The latter in turn yields : \y\ ~ (1/18 to 1/23.6), with y having the same sign as 
e (see Eq. (|3)). This solution for y obtains only by assuming that y is 0(1/10) rather than 
0(1). Combining now with the mixing in the /i-r sector determined above (see Eq. (P^), 
one can then determine the z/^-z/,- oscillation angle. The two predictions of the model for the 
neutrino-system are then : 

m(z/^) ^ (1/10- 1/30) eV (29) 



C ^ O^r-O^r ^ ^0.437 + ^-^ J . (30) 

Thus, sin^ 261^1^ = (0.96, 0.91, 0.86, 0.83, 0.81) (31) 
for m^Jm^^ = (1/10, 1/15, 1/20, 1/25, 1/30) . (32) 

Both of these predictions are extremely successful. 

Note the interesting point that the MSW solution, together with the requirement that \y\ 
should have a natural hierarchical value (as mentioned above) , lead to y having the same sign 
as e; that (it turns out) implies that the two contributions in Eq.(|30|) must add rather than 



subtract, leading to an almost maximal oscillation angle [14]. The other factor contributing 
to the enhancement of 6*°^^^ is, of course, also the asymmetry-ratio which increases |^^^| 
from 0.25 to 0.437 (see Eq. (0)). We see that one can derive rather plausibly a large i'fj,-i'T 
oscillation angle sin^ 2^°^^^ > 0.8, together with an understanding of hierarchical masses 
and mixings of the quarks and the charged leptons, while maintaining a large hierarchy in 



This estimate for Mr is retained even if one allows for i^^-iyr mixing (see Ref. IT 
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the seesaw derived neutrino masses [mi^^/m,^^ = 1/10-1/30), all within a unified framework 
including both quarks and leptons. In the example exhibited here, the mixing angles for 
the mass eigenstates of neither the neutrinos nor the charged leptons are really large, in 
that ef^^ ~ 0.437 ~ 23° and 9''^^ ~ (0.18-0.31) ^ (10-18)°, yet the oscillation angle obtained 
by combining the two is near-maximal. This contrasts with most works in the literature 
in which a large oscillation angle is obtained either entirely from the neutrino sector (with 
nearly degenerate neutrinos) or almost entirely from the charged lepton sector. 

While Mji ^ (5-15) x 10^^ GeV and y ^ —1/20 are better determined, the parameters x 
and z can not be obtained reliably at present because very little is known about observables 
involving i/g- Taking, for concreteness, m^^ ^ (10~^-10~^ (1 to few)) eV and 6°^'^ ~ d^^ — O'^^ ^ 
10^^ ±0.03 as inputs, we obtain : z ~ (1-5) x 10^^ and x ~ (1 to few)(10^^-10~^), in accord 
with the guidelines of 1^1 ~ ||/|/10 and |a;| ~ z"^. This in turn yields : 6°^^^ ~ — 6'^^ ~ 
0.06±0.015. Note that the mass of m^^ ~ 3 x 10"'^ eV, that follows from a natural hierarchical 
value for y ~ —(1/20), and Ben as above, go well with the small angle MSW explanation^ of 
the solar neutrinos puzzle. 

It is worthnoting that although the superheavy Majorana masses of the RH neutrinos 
cannot be observed directly, they can be of cosmological significance. The pattern given 
above and the arguments given in Sec.|T| and in this section suggests that M(i/]j) ^ (5-15) x 
10^^ GeV, M(z/^) ^ (1-4) x 10^^ Q^y (for x ^ 1/20); and M(^'^) ~ (1/2-10) x 10^ GeV (for 
X ~ (l/2-10)10~^ > z^). A mass of z/^ ~ 10^ GeV is of the right magnitude for producing 
z/|. following reheating and inducing lepton asymmetry in i/Jj decay into + u}^, that is 
subsequently converted into baryon asymmetry by the electroweak sphalerons |T^JT7[] . 

In summary, we have proposed an economical and predictive pattern for the Dirac mass 
matrices, within the SO(10)/G(224)-framework, which is remarkably successful in describing 
the observed masses and mixings of all the quarks and charged leptons. It leads to five 
predictions for just the quark- system, all of which agree with observation to within 10%. 
The same pattern, supplemented with a similar structure for the Majorana mass matrix, 
accounts for both the large z/^-z/^ oscillation angle and a mass of z/^ ~ 1/20 eV, suggested 
by the SuperK data. Given this degree of success, it makes good sense to study proton 
decay concretely within this SO(10)/G(224)-framework. The results of this study JI^ are 
presented in the next section. 

Before turning to proton decay, it is worth noting that much of our discussion of fermion 
masses and mixings, including those of the neutrinos, is essentially unaltered if we go to the 
limit e' — of Eq. (28). This limit clearly involves: 

= 0, Be ^ a/ nid/ms , m^^ = , 6^^ = 6'^^ = . 



\Vub\ ^ \l ^^—^ V^d/rrib (ms/mb) - (2.1)(0.039)(0.023) ~ 0.0019 (33) 



^Although the small angle MSW solution appears to be more generic within the approach outlined 
above, we have found that the large angle solution can still plausibly emerge in a limited region of 
parameter space, without affecting our results on fermion masses. 
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All other predictions remain unaltered. Now, among the observed quantities in the list above, 
6c — \fmsjrns is a good result. Considering that mu/mt ~ 10~^, = is also a pretty 
good result. There are of course plausible small corrections which could arise through Planck 
scale physics; these could induce a small value for through the (l,l)-entry 5 ~ 10~^. For 
considerations of proton decay, it is worth distinguishing between these two extreme variants 
which we will refer to clS CcLSGS I and II respectively. 

Case I :e' ^ 2 X 10"^ ,5 = 

Case II -5 ^ 10"^ , e' = . (34) 

It is worth noting that the observed value of \Vub\ ~ 0.003 favors a non-zero value of e'(^ 
(1 — 2) X 10~^). Thus, in reality, e' may not be zero, but it may lie in between the two 
extreme values listed above. In this case, the predicted proton lifetime for the standard 
d = 5 operators would be intermediate between those for the two cases, presented in Sec.|VI[ 



VI. EXPECTATIONS FOR PROTON DECAY IN SUPERSYMMETRIC UNIFIED 

THEORIES 



A. Preliminaries 



Turning to the main purpose of this talk, I present now the reason why the unification 
framework based on SUSY SO(IO) or G(224), together with the understanding of fermion 
masses and mixings discussed above, strongly suggest that proton decay should be imminent. 

Recall that supersymmetric unified theories (GUTs) introduce two new features to proton 
decay : (i) First, by raising Mx to a higher value of about 2 x 10^® GeV (contrast with the 
non-supersymmetric case of nearly 3 x lO^'' GeV), they strongly suppress the gauge-boson- 
mediated d = 6 proton decay operators, for which e~^7r^ would have been the dominant mode 
(for this case, one typically obtains : T^^{p e+7r*')|rf=6 ~ io35-3±i-5 j^s). (ii) Second, they 
generate d = 5 proton decay operators of the form QiQjQkQi/M in the superpotential, 
through the exchange of color triplet Higginos, which are the GUT partners of the standard 
Higgs(ino) doublets, such as those in the 5-1-5 of SU(5) or the 10 of SO(IO). Assuming that a 
suitable doublet-triplet splitting mechanism provides heavy GUT-scale masses to these color 
triplets and at the same time light masses to the doublets, these "standard" d = 5 operators, 
suppressed by just one power of the heavy mass and the small Yukawa couplings, are found 



to provide the dominant mechanism for proton decay in supersymmetric GUT [^^ -|47[|. 

Now, owing to (a) Bose symmetry of the superfields in QQQL/M, (b) color antisymmetry, 
and especially (c) the hierarchical Yukawa couplings of the Higgs doublets, it turns out that 
these standard d = 5 operators lead to dominant VK^ and comparable T/tt"*" modes, but in 
all cases to highly suppressed e~^n^, e~^K^ and even ^^K^ modes. For instance, for minimal 
SUSY SU(5), one obtains (with tan/5 < 20, say) : 

[r(/i+if°)/r(i7K+)]f,^(') ~ K/(m, siYi^e,)fR ^ lo-^, (35) 

where i? ~ 0.1 is the ratio of the relevant |matrix element px (phase space), for the two 
modes. 
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It was recently pointed out that in SUSY unified theories based on SO(IO) or G(224), 
which assign heavy Majorana masses to the RH neutrinos, there exists a new set of color 
triplets and thereby very likely a new source of c? = 5 proton decay operators |T^. For 
instance, in the context of the minimal set of Higgs multiplets^^ {45/;-, 16h, 16^ and 10^} 
(see Sec.[V]), these new d = 5 operators arise by combining three effective couplings intro- 
duced before :- i.e., (a) the couplings /jjl6jl6jl6//16H/M (see Eq.(^) that are required 
to assign Majorana masses to the RH neutrinos, (b) the couplings gijlGilGjlGn^GH /M , 
which are needed to generate non-trivial CKM mixings (see Eq.(|15D), and (c) the mass term 
MiqIGh^Qh- For the fij couplings, there are two possible SO(10)-contractions (leading to 
a 45 or a 1) for the pair 16jl6//, both of which contribute to the Majorana masses of the 
RH neutrinos, but only the non-singlet contraction (leading to 45), would contribute to d=5 
proton decay operator. In the presence of non-perturbative quantum gravity, one would in 
general expect the two contractions to have comparable strength. Furthermore, the couplings 
of 45 's lying in the string-tower or possibly below the string-scale, and likewise of singlets, to 
the 16j- IG^-pair, would respectively generate the two contractions. It thus seems most likely 
that both contractions are present, having comparable strength. Allowing for a difference 
between the relevant projection factors for masses versus proton decay, and also for the 
fact that both contractions contribute to the former, but only the non-singlet one (i.e. 45) 
to the latter, we would set the relevant fij coupling for proton decay to be {fij)p = {fij)u ■ K, 
where {fij)u defined in Sec.|V| directly yields - masses (see Eq.(|D); and K is a relative 
factor of order unity. As a plausible range, we will take 1/3 to 2 (say). In the presence 
of the non-singlet contraction, the color-triplet Higginos in 16h and IQh of mass Mig can be 
exchanged between qiqj and g^grpairs (correspondingly, for G(224), the color triplets would 
arise from (1, 2, 4:)h and (1, 2, 4:)h). This exchange generates a new set of d = 5 operators in 
the superpotential of the form 

Wne^ cc if,,), qmK {16,lQj) (IQklQi) {TQh) {IQh)/M^ X (l/Mig), (36) 

which induce proton decay. Note that these operators depend, through the couplings fij 
and Qki, both on the Majorana and on the Dirac masses of the respective fermions. This is 
why within SUSY SO(IO) or G(224), proton decay gets intimately linked to the masses and 
mixings of all fermions, including neutrinos. 

B. Framework for Calculating Proton Decay Rate 

To establish notations, consider the case of minimal SUSY SU(5) and, as an example, the 
process cd su^, which induces p —>■ V^K~^. Let the strength of the corresponding d = 5 
operator, multiplied by the product of the CKM mixing elements entering into wino-exchange 
vertices, (which in this case is sin 6^(7 cos 6^(7) be denoted by A. Thus (putting cos6'c = 1), 
one obtains: 



The origin of the new d = 5 operators in the context of other Higgs multiplets, in particular in 
the cases where 126h and 126// are used to break B-L, has been discussed in Ref. p!9| ]. 
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A,^{SU{5)) = {KMMhc) sin^e ^ {m,m, sin^ 0^/^') (tan/^/M^J 

(1.9 X 10"®) (tan/?/MHe) ^ (2 x lO^^^GeV-^) (tan/3/2) (2 x If}^"^ GeV / Mhc, 



(37) 



where tan/3 = Vujvd-, and we have put f„ = 174 GeV and the fermion masses extrapolated 
to the unification-scale - i.e. mc ~ 300 MeV and ~ 40 MeV. The amplitude for the 
associated four-fermion process dus — * z7^ is given by: 



A^{dus ^ V,:) = X (2/) 



(3^ 



where / is the loop-factor associated with wino-dressing. Assuming rriyj -C rriq ~ m^-, one 
gets: / ~ (miu/m?) (a2/47r). Using the amplitude for ((i'u)(sz/^), as in Eq. (|38|), ( 



/i or T) 



one then obtains 14 



r-\p VrK+) ^ (0.6 X 10^^)yrs x 



0.67 

'a7 



0.014 GeV^'^ 



H 



(1/6) 



_(m^^/mg-) 



1.2 TeV, 



2 X 10-^4 GeV 



T-1 



n 2 



(39) 



Here /?// denotes the hadronic matrix element defined by PnUiik) = ea/3'y{0\{d1u^)ul\p, k) . 
While the range Ph = (0.003-0.03) GeV^ has been used in the past |4^, given that one 
lattice calculation yields (3h = (5.6±0.5) x 10~^ GeV^ [Q, and a recent improved calculation 
yields ~ 0.014GeV^ (whose systematic errors that may arise from scaling violations 
and quenching are hard to estimate |^), we will take as a conservative, but plausible, 
range for Ph to be given by (0.014GeV'^)(l/2 — 2). [Compare this with the range for = 
(0.006GeV3)(l/2 - 2) as used in Ref. 11]. Here, As ^ 0.67 stands for the short distance 
renormalization factor of the d = 5 operator. Note that the familiar factors that appear in 
the expression for proton lifetime - i.e., Mhci (1 + 2/^) representing the interference between 



the t and c contributions, and tan/3 (see e.g. Ref. |^ and discussion in the Appendix 
of Ref. H^) - are all effectively contained in A(z7). In Ref. |14|, guided by the demand of 
naturalness (i.e. absence of excessive fine tunning) in obtaining the Higgs boson mass, squark 
masses were assumed to lie in the range of 1 TeV(l/y2 - 72), so that m^- < 1.4TeV. Recent 
work, based on the notion of focus point supersymmetry however suggests that squarks may 
be considerably heavier without conflicting with the demands of naturalness ||50|. In the 
interest of obtaining a conservative upper limit on proton lifetime, we will therefore allow 
squark masses to be as heavy as about 2.5 TeV and as light as perhaps 600 GeV. 



We remark that if the recently reported (g-2) - anomaly for the muon |51| is attributed to 
supersymmetry |52|, one would need to have extremely light s-fermions (i.e. m^- ~ 200 — 400 GeV 
(say) and correspondingly (for promising mechanisms of SUSY-breaking) rriq < 300 — 600 TeV 
(say)), and simultaneously large or very large tan /3(~ 25 — 50). However, not worring about grand 
unification, such light s-fermions, together with large or very large tan/3 would typically be in 
gross conflict with the limits on the edm's of the neutron and the electron, unless on can explain 
naturally the occurence of minuscule phases (< 1/300 to 1/1000) and/or large cancellation. Thus, 
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Allowing for plausible and rather generous uncertainties in the matrix element and the 
spectrum we take: 

/3h = (0.014 GeV^) (1/2-2) 



m^/rrig = 1/6(1/2-2), and niq ^ m^- 1.2 TeV (1/2 -2) . (40) 
Using Eqs.(|39l-|0D, we get: 



T-\p V^K+) ^ (0.6 X 10^Vrs)[2x 10"^^ GeV-7A(l7^) ]^ {64- 1/64 } . (41) 

Note that the curly bracket would acquire its upper-end value of 64, which would serve to- 
wards maximizing proton lifetime, only provided all the uncertainties in Eq. (p]) are stretched 
to the extreme so that (3h = 0.007 GeV^, m^r/mq ^ 1/12 and rriq ^ 2.4 TeV. This relation, 
as well as Eq. ( pOf ) are general, depending only on Aiyi) and on the range of parameters 
given in Eq. (140|) . They can thus be used for both SU(5) and SO(IO). 

The experimental lower limit on the inverse rate for the uK~^ modes is given by 

[J2 r(p - T^eK+)]-j^, > 1.6 X lO^^'yrs. (42) 

e 

Allowing for all the uncertainties to stretch in the same direction (in this case, the curly 
bracket = 64), and assuming that just one neutrino flavor (e.g. z/^ for SU(5)) dominates, the 
observed hmit (Eq.(^)) provides an upper bound on the amphtudeQ: 

A{Ve) < 1 X 10-24 GeV"^ (43) 

which holds for both SU(5) and SO(IO). Recent theoretical analyses based on LEP-limit on 
Higgs mass (> 114 GeV), together with certain assumptions about MSSM parameters (as 



if the {g — 2)^ - anomaly turns out to be real, it may quite possibly need a non-supersymmetric 
explanation, in accord with the edm-constraints which ordinarily seem to suggest that squarks are 
(at least) moderately heavy {rriq > 0.6— 1 TeV, say), and tan/3 is not too large (< 3 to 10, say). We 
mention in passing that the extra vector - like matter - specially a 16 + 16 of SO(IO) - as proposed 



in the so-called extended supersymmetric standard model (ESSM) |^,53|, with the heavy lepton 
mass being of order (150-200) hundred GeV, can provide such an explanation |54]. Motivations for 
the case of ESSM, based on the need for (a) removing the mismatch between MSSM and string 
unification scales, and (b) dilaton-stabilization, have been noted in Ref. |^^. Since ESSM is an 
interesting and viable variant of MSSM, and would have important implications for proton decay, 
we will present the results for expected proton decay rates for the cases of both MSSM and ESSM 
in the discussion to follow. 

^^If there are sub-dominant T'iK'^ modes with branching ratio R, the right side of Eq. (^) should 
be divided by \/l + R. 
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in CMSSM) and/or constaint from muon g-2 anomaly |5T[] suggest that tan/3 > 3 to 5 |[56|| . 
In the interest of getting a conservative upper limit on proton lifetime, we will therefore use, 
as a conservative lower limit, tan/3 > 3. We will however exhibit relevant results often as 
a function of tan (3 and exhibit proton lifetimes corresponding to higher values of tan/5 as 
well. For minimal SU(5), using Eq.(0) and, conservatively tan/3 > 3, one obtains a lower 
limit on Mfjc given by: 

Mhc > 5.5 X lO^^GeV (SU(5)) (44) 

At the same time, higher values of Mhc > 3 x 10^^ GeV do not go very well with gauge 
coupling unification. Thus we already see a conflict, in the case of minimal SUSY SU(5), 
between the experimental limit on proton lifetime on the one hand, and coupling unification 
and constraint on tan (3 on the other hand. To see this conflict another way, if we keep Mhc ^ 
3 X 10^^ GeV (for the sake of coupling unification) we obtain from Eq.(^): A(SU(5)) > 
1.9 X IQ-^^ GeV"^(tan/3/3). Using Eq. (|T]), this in turn implies that 

r-^(j9 ^ z7A^+) < 0.6 X lO^^yrs x (3/tan/5)^ (SU(5)) (45) 

For tan/3 > 3, a lifetime of 0.7 x 10^^ years is thus a conservative upper limit. In practice, it 
is unlikely that all the uncertainties, including these in Mhc and tan/3, would stretch in the 
same direction to nearly extreme values so as to prolong proton lifetime. A more reasonable 
upper limit, for minimal SU(5), thus seems to be: T~^{p — > VK~^){SU{5)) < (0.3) x 10^^ yrs. 
Given the experimental lower limit (Eq.(|42|)), we see that minimal SUSY SU(5) is already 
excluded (or strongly disfavored) by proton decay-searches. We have of course noted in 
Sec.|V] that SUSY SU(5) does not go well with neutrino oscillations observed at SuperK. 

Now, to discuss proton decay in the context of supersymmetric SO(IO), it is necessary 
to discuss first the mechanism for doublet-triplet splitting. Details of this discussion may be 



found in Ref. |14]. A synopsis is presented in the Appendix. 



C. Proton Decay in Supersymmetric SO(IO) 

The calculation of the amplitudes Agtd and Anew for the standard and the new operators 
for the SO(IO) model, are given in detail in Ref. Here, I will present only the results. It 
is found that the four amplitudes Agtdi'J^T-K'^), Astdi'^iJ.K'^), Anew^prK^) and AnewiJ'^lK^) 
are in fact very comparable to each other, within about a factor of two to five, either way. 
Since there is no reason to expect a near cancellation between the standard and the new 
operators, especially for both V^-K^ and V^K^ modes, we expect the net amplitude (standard 
-|- new) to be in the range exhibited by either one. Following Ref. [Q, I therefore present 
the contributions from the standard and the new operators separately. 

One important consequence of the doublet-triplet splitting mechanism for SO(IO) out- 
lined briefly in the appendix and in more detail in Ref. is that the standard d=5 proton 
decay operators become inversely proportional to Mf,ff = [A < A^h >]^/^io' ~ M'j^/MiQi, 
rather than to Mh^- Here, Mio' represents the mass of 10'^^^, that enters into the D-T split- 
ting mechanism through effective coupling A10h45h10^ in the superpotential (see Appendix, 
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Eq.([M|)). As noted in Ref. [|I^], Miqi can be naturally suppressed (due to flavor symme- 
tries) compared to Mx, and thus Me// correspondingly larger than Mx by even one to 
three orders of magnitude. It should be stressed that Me// does not represent the physical 
masses of the color triplets or of the other particles in the theory. It is simply a parameter 
of order Mjr/MiQi. Thus larger values of M^jf, close to or even exceeding the Plank scale, 
do not in any way imply large corrections from quantum gravity. Now accompanying the 
suppression due to Me//, the standard proton decay amplitudes for SO(IO) possess an in- 
trinsic enhancement as well, compared to those for SU(5), owing primarily due to differences 



in their Yukawa couplings for the up sector (see Appendix C of Ref. |jTj])- As a result of 
this enhancement, combined with the suppression due to higher values of Mg//, a typical 
standard d = 5 amphtude for SO(IO) is given by (see Appendix C of Ref. fl^ 



M^.K^fsT'^ ^ (/^33/Me//)(2 X 10-^), 

which should be compared with A{iJf,K+ff/^ ^ (1.9 x 10-®)(tan/5/MjyJ (see Eq.{^). 
Note, taking /133 ^ 1/4, the ratio of a typical SO(IO) over SU(5) amplitude is given by 
(Mii/^/Me//)(88)(3/ tan/3). Thus the enhancement by a factor of about 88 (for tan/5 = 3), 
of the SO(IO) compared to the SU(5) amplitude, is compensated in part by the suppression 
that arises from Me// being larger than M^^- 

In addition, note that in contrast to the case of SU(5), the SO(IO) amplitude does not 
depend explicitly on tan j3. The reason is this: if the fermions acquire masses only through 
the IOh in SO(IO), as is well known, the up and down quark Yukawa couplings will be 
equal. By itself, it would lead to a large value of tan P = rrit/mb ^ 60 and thereby to a large 
enhancement in proton decay amplitude. Furthermore, it would also lead to the bad relations: 
rric/rns = mt/nif, and Vckm = 1- However, in the presence of additional Higgs multiplets, 
in particular with the mixing of {16H)d with IOh (see Appendix and Sec.|V|), (a) tan/3 can 
get lowered to values like 3-20, (b) fermion masses get contributions from both < 16^ >d 
and < IOh >, which correct all the bad relations stated above, and simultaneously (c) the 
explicit dependence of A on tan (3 disappears. It reappears, however, through restriction on 
threshold corrections, discussed below. 

Although Me// can far exceed Mx, it still gets bounded from above by demanding that 
coupling unification, as observed 0, should emerge as a natural prediction of the theory 
as opposed to being fortuitous. That in turn requires that there be no large (unpredicted) 
cancellation between GUT-scale threshold corrections to the gauge couplings that arise from 
splittings within different multiplets as well as from Plank scale physics. Following this point 
of view, we have argued (see Appendix) that the net "other" threshold corrections to a^lniz) 



13 Pqj. instance, in the absence of GUT-scale threshold corrections, the MSSM value of 
Oi3{mz)MSSM, assuming coupling unification, is given by ck3("T'z)m5SM ~ 0.125 — 0.13 which 
is about 5-8% higher than the observed value: Oi3{'>nz)°MssM ~ 0-118 — 0.003 [^]. We demand 
that this discrepancy should be accounted for accurately by a net negative contribution from D-T 
splitting and from "other" threshold corrections (see Appendix, Eq.(|A4|)), without involving large 
cancellations. That in fact does happen for the minimal Higgs system (45, 16, 16) [see Ref. |14]]. 
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arising from the Higgs (in our case 45^, 16// and IQh) and the gauge multiplets should be 
negative, but conservatively and quite plausibly no more than about 10%. This in turn 
restricts how big can be the threshold corrections to 0:3 (m^) that arise from (D-T) splitting 
(which is positive). Since the latter is proportional to ln(Me// cos7/Mx) (see Appendix), we 
thus obtain an upper limit on Mg// cos 7. For the simplest model of D-T splitting presented 
in Ref. [|14| and in the Appendix (Eq. ([Al|) ), one obtains: COS7 ~ {ta.n j3) / {rrit / rrih) . An upper 
limit on Me// cos 7 thus provides an upper limit on Me// which is inversely proportional to 
tan (3. In short, our demand of natural coupling unification, together with the simplest model 
of D-T splitting, introduces an implicit dependence on tan f3 into the lower limit of the SO(IO) 
- amplitude - i.e. A{SO{10)) oc 1/Me// > (a quantity) oc tan/3. These considerations are 
reflected in the results given below. 

Assuming tan/3 > 3 and accurate coupling unification (as described above), one obtains 
for the case of MSSM, a conservative upper limit on Me// < 2.7 x 10^^ GeV (3/ tan/?) (see 
Appendix and Ref. |]I1|). Using this upper limit, we obtain a lower limit for the standard 
proton decay amplitude given by 



std 



> 



(7.8 X 10-24 GeV-^) (1/6 - 1/4) easel 
(3.3 X 10-24 GeV-^) (1/6 - 1/2) case II 



/ SO(10)/MSSM, with 
\^ tan /3 > 3 



(46) 



Substituting into Eq.(^) and adding the contribution from the second competing mode 



Ui^K^ ^ with a typical branching ratio R ~ 0.3, we obtain 



T-\vK+)std < 



(0.7 X lO^'Vrs.) (1.6-0.7) 
(1.5 X lO^Vrs.) (4-0.44) 



{64-1/64} 



SO(10)/MSSM, with\ 
tan/3 > 3 J' 



(47) 



The upper and lower entries in Eqs.(^6D and (^7]) correspond to the cases I and II of the 
fermion mass-matrix with the extreme values of e' - i.e. e' = 2 x 10-"^ and e' = - respectively, 
(see Eq.(0)). The uncertainty shown inside the square brackets correspond to that in the 
relative phases of the different contributions. The uncertainty of {64 to 1/32} arises from 
that in /?//, (m^y/m^) and rriq (see Eq.(^OD). Thus we find that for MSSM embedded in 
SO (10), for the two extreme values of e' (cases I and II) as mentioned above, the inverse 
partial proton decay rate should satisfy: 



_ ni+2.0 

0.7 X 10^ -i 7yrs. 

yrs. 



1.3 X lO^^-i'se 



< 



0.7 X 10^3 yrs. 
3.7 X 10^3 yrs. 



SO(10)/MSSM, with 
tan /? > 3 



The central value of the upper limit in Eq.(^S]) corresponds to taking the upper limit on 
^eff < 2.7 X 10^^ GeV, which is obtained by restricting threshold corrections as described 
above (and in the Appendix) and by setting (conservatively) tan/3 > 3. The uncertainties of 
matrix element, spectrum and choice of phases are reflected in the exponents. The uncertainty 
in the most sensitive entry of the fermion mass matrix - i.e. e' - is fully incorporated (as 
regards obtaining an upper limit on the lifetime) by going from case I (with e' = 2 x 10-^) 
to case II (e' = 0). Note that this increases the lifetime by almost a factor of six. Any 
non-vanishing intermediate value of e' would only shorten the lifetime compared to case II. 
In this sense, the larger of the two upper limits quoted above is rather conservative. We see 
that the predicted upper limit for case I of MSSM (with the extreme value of e' = 2 x 10^^) 
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is already in conflict with the empirical lower limit (Eq.(^)) while that for case II i.e. e' = 
(with all the uncertainties stretched as mentioned above) is only about two times higher 
than the empirical limit. 

Thus the case of MSSM embedded in SO(IO) is already tightly constrained, to the point 
of being rather disfavored, by the limit on proton lifetime in that all the parameters need 
to lie near their "extreme" ends so that it may be compatible with the empirical limit (see 
also results for other choices of parameters listed in Table 1). The constraint is of course 
augmented especially by our requirement of natural coupling unification which prohibits 
accidental large cancellation between different threshold corrections (see Appendix); and 
it will be even more severe, especially within the simplest mechanism of D-T splitting (as 
discussed in the Appendix), if tan P turns out to be larger than 5 (say). On the positive side, 
improvement in the current limit by a factor of even 2 to 3 ought to reveal proton decay, 
otherwise the case of MSSM embedded in SO (10), would be clearly excluded. 



D. The case of ESSM 

Before discussing the contribution of the new d = 5 operators to proton decay, an in- 
teresting possibility, mentioned in the introduction, that would be especially relevant in the 
context of proton decay, if tan P is large, is worth noting. This is the case of the extended 
super symmetric standard model (ESSM), which introduces an extra pair of vector-like fam- 
ilies (16 + I6 of SO(IO)), at the TeV scale pOlJSl. Adding such complete SO(10)-muhiplets 



would of course preserve coupling unification. From the point of view of adding extra fami- 
lies, ESSM seems to be the minimal and also the maximal extension of the MSSM, that is 
allowed in that it is compatible with (a) neutrino-counting, (b) precision electroweak tests, 
as well as (c) a semi-perturbative as opposed to non-perturbative gauge coupling unification 
20| , |53[] . Q The existence of two extra vector-like families can of course be tested at the LHC. 

Theoretical motivations for the case of ESSM arise because, (a) it raises Ounif to a semi- 
perturbative value of 0.25 to 0.3, and therefore has a better chance to achieve dilaton- 
stabilization than the case of MSSM, for which Ounif is rather weak (only 0.04); and (b) owing 
to increased two-loop effects |^,S7|, it raises the unification scale Mx to (1/2 — 2) x 10^''GeV 



and thereby considerably reduces the problem of a mismatch between the MSSM and the 
string unification scales (see Sec.|ITT|). A third feature relevant to proton decay is the following. 
In the absence of unification-scale threshold and Planck-scale effects, the ESSM value of 
O-zijnz) obtained by assuming gauge coupling unification, which we denote by a3(?^z)EssM 
is lowered to about 0.112 — 0.118 compared to a3("^z)MSSM ~ 0.125 — 0.13. 



As explained in the appendix, the net result of these two effects - i.e. a raising of Mx and 
a lowering of a3('TT'z)ESSM " that for ESSM embedded in SO(IO), tan/3 can span a wide 
range from 3 to even 30, and simultaneously the value or the upper limit on Me// can range 
from (60 to 6) x lO^^GeV, in full accord with our criterion for accurate coupling unification 



For instance, addition of two pairs of vector-like families at the TeV-scale, to the three chiral 
families, would cause gauge couplings to become non-perturbative below the unification scale. 
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discussed above. 

Thus, in contrast to MSSM, ESSM allows for larger values of tan (3 (like 20 or 30), without 
needing large threshold corrections, and simultaneously without conflicting with the limit 
on proton lifetime. 

To be specific, consider first the case of a moderately large tan/? = 20 (say), for which 
one obtains Me// ~ 9 x 10^^ GeV, with the "other" threshold correction — being about 
5% (see Appendix for definition). In this case, one obtains: 



Istd 



'1.6-0.7) 
(10-1) 



{64 - 1/64} (7 X 103^ yrs) ( ^''^'Zfl'to ) ' ^''^ 



As before, the upper and lower entries correspond to cases I (e' = 2 x lO"'') and II (e' = 0) 
of the fermion mass- matrix (see Eq.(|^)). The uncertainty in the upper and lower entries 
in the square bracket of Eq . (^91) corresponds to that in the relative phases of the different 
contributions for the cases I and II respectively, while the factor {64-1/64} corresponds to 
uncertainties in the SUSY spectrum and the matrix element (see Eq.^). 

We see that by allowing for an uncertainty of a factor of (30 — 100) jointly from the two 
brackets for Case I (and (13 — 44) for Case II), proton lifetime arising from the standard 
operators would be expected to lie in the range of (2.2 — 7.5) x 10^'^ yrs, for the case of ESSM 
embedded in SO(IO), with tan/5 = 20. Such a range is compatible with present limits, but 
accessible to searches in the near future. 

The other most important feature of ESSM is that, by allowing for larger values of 
Me//, especially for smaller values of tan/? 3 to 10 (say), the contribution of the standard 
operators by itself can be perfectly consistent with present limit on proton lifetime even for 
almost central or "median" values of the parameters pertaining to the SUSY spectrum, the 
relevant matrix element, e' and the phase- dependent factor. 

For instance, for ESSM, one obtains Me// f» (4.5 x GeV) {A /tan (3), with the "other" 
threshold correction - 6'^ being about 5% (see Appendix and Eq. ( |A6| ) ) . Now, combining 
cases I (e' = 2 X 10^^) and II (e' = 0), we see that the square bracket in Eq.(|49|) which we 
will denote by [S], varies from 0.7 to 10, depending upon the relative phases of the different 
contributions and the values of e'. Thus as a "median" value, we will take [^jmed ~ 2 to 6. 
The curly bracket {64-1/64}, to be denoted by {C}, represents the uncertainty in the SUSY 
spectrum and the matrix element (see Eq. (^OD ) . Again as a "nearly central" or "median" 
value, we will take {C}med ~ 1/6 to 6. Setting Me// as above we obtain 

r-i(i^ir+);7'^-"" ^ [5]„ed{C}w(1.8 X 10^3 yrs)(4/tan/3)2(SO(10)/ESSM). (50) 

Choosing a few sample values of the effective parameters [S] and {C}, with low values of 
tan/3 = 4 to 10, the corresponding values of T~^{h'K~^), following from Eq . (|50D , are listed 
below in Table 1. 

Note that ignoring contributions from the new d=5 operators for a moment the entries 
in Table 1 represent a very plausible range of values for the proton lifetime, for the case of 
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As I will discuss in the next section, we of course expect the new d=5 operators to be important 
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ESSM embedded in SO(IO), with tan (3 ^ 3 to 10 (say), rather than upper limits for the 
same. This is because they are obtained for "nearly central" or "median" values of 



TABLE 1. PROTON LIFETIME, BASED ON CONTRIBUTIONS FROM 
ONLY THE STANDARD OPERATORS FOR THE CASE OF ESSM 
EMBEDDED IN SO(IO), WITH PARAMETERS BEING IN THE 

"MEDIAN" RANGE 



tan l3 = A 
[S]=2.7 
{C}=l/2 to 2 


tan /? = 4 
[S]=6 
{C}=l/6 to 1 


tan /5 = 10 

[S]=5.4 
{C}=1 to 6 


tan/? = 10 

[S]=6 
{C}=1 to 4 


r-'iuK+)rssM^ 
(2.5 to 10) X 10=^^ yrs 


(1.8 to 11) X 10^^ yrs 


(1.6 to 10) X 10^3 yrs 


(1.8 to 7.3) X 10=^^ yrs 



the parameters represented by the values of [S]~ 2 to 6 and {C}~ 1/6 to 6, as discussed 
above. For instance, consider the cases {C}=1 and {C}=l/6 respectively, both of which (as 
may be inferred from the table) can quite plausibly yield proton lifetimes in the range of 

= 0.014GeV^ (the central value of 
I)), while that of {C}=l/6 would 
i 600GeV and m^^/rrig 1/5. In 



10 to 10 yrs. Now {C}=1 corresponds, e.g., to (3h 
1.2 TeV and my^/rriq = 1/6 (see Eq.([ 
0.014GeV^ with rrin 



Ref. H 



rUg 



correspond, for example, to Ph 
short, for the case of ESSM, with low values of tan P ^ 3 to 10 (say), squark masses can be 
well below 1 TeV, without conflicting with present limit on proton lifetime. This feature is 
not permissible within MSSM embedded in SO(IO). 

Thus, confining for a moment to the standard operators only, if ESSM represents low- 
energy physics, and if tan/3 is rather small (3 to 10, say), we do not have to stretch at all 
the uncertainties in the SUSY spectrum and the matrix elements to their extreme values (in 
contrast to the case of MSSM) in order to understand why proton decay has not been seen 
as yet, and still can be optimistic that it ought to be discovered in the near future, with a 
lifetime < 10^^ years. The results for a wider variation of the parameters are listed in Table 
2, where contributions of the new d=5 operators are also shown. 

It should also be remarked that if in the unlikely event, all the parameters (i.e. (5hi 
{my^/rriq), rriq and the phase-dependent factor) happen to be closer to their extreme values 
so as to extend proton lifetime, and if tan (3 is small (~ 3 to 10, say) and at the same time the 
value of Me// is close to its allowed upper limit (see Appendix), the standard d=5 operators 
by themselves would tend to yield proton lifetimes exceeding even (1/3 to l)xlO^^ years for 
the case of ESSM, (see Eq.(^) and Table 2). In this case (with the parameters having nearly 
extreme values), however, as I will discuss shortly, the contribution of the new d=5 operators 
related to neutrino masses (see Eq.(^)), would dominate and quite naturally yield lifetimes 
bounded above in the range of (1 — 10) x 10^^ years (see Sec. |VI E| and Table 2). Thus in the 



and significantly infuence proton lifetime (see e.g. Table 2). Entries in Table 1 could still represent 
the actual expected values of proton lifetimes, however, if the parameter K defined in VI A| (also 
see VI E|) happens to be unexpectedly small (^ 1). 
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presence of the new operators, the range of (10^^ — 10^^) years for proton lifetime is not only 
very plausible but it also provides a conservative upper limit, for the case of ESSM embedded 
m SO(IO). 



E. Contribution from the new d=5 operators 



As mentioned in Sec. [VI A] , for supersymmetric G(224)/SO(10), there very hkely exists 



a new set of d=5 operators, related to neutrino masses, which can induce proton decay 
(see, Eq.(|4^)). The decay amplitude for these operators for the leading mode (which in this 
case is v^K^) becomes proportional to the quantity P = {{f33)u{16H) /M}h33K/ {Miq tan 'y), 
where (/33)jy and /133 are the effective couplings defined in Eqs.(0) and (|T5|) respectively, and 
M16 and tan 7 are defined in the Appendix. The factor K, defined by (/33)p = {f2,z)yK, is 
expected to be of order unity (see SecjV^A for the origin of K). As a plausible range, we 
take K ^ 1/3 to 2. Using Migtan7 = X'{16h) (see Appendix), and /133 ~ 1/2 (given by top 
mass), one gets: P ~ {{f3s)iy/M){l/2X')K. Here M denotes the string or the Planck scale 
(see Sec.|I3 and footnote 2); thus M pa (1/2 — 1) x lO^'^GeV; and A' is a quartic coupling 



defined in the appendix. Validity of perturbative calculation suggests that A' should not 
much exceed unity, while other considerations suggest that A' should not be much less than 
unity either (see Ref. []T1|, Sec.6E). Thus, a plausible range for A' is given by A' ~ (1/2 — y^). 



(Note it is only the upper limit on A' that is relevant to obtaining an upper limit on proton 
lifetime). Finally, from consideration of Ur mass, we have {f33)u ~ 1 (see Sec.|Vp. We thus 
obtain: P ^ (5x 10~^^GeV~^)(l/-\/2 to 4)K. Incorporating a further uncertainty by a factor 
of (1/2 to 2) that arises due to choice of the relative phases of the different contributions 
(see Ref. ^^), the effective amplitude for the new operator is given by 

A{Uf,K+)^ew ~ (1.5 X 10-24GeV-^)(l/2V2 to 8)K (51) 

Note that this new contribution is independent of M^ff] thus it is the same for ESSM as 
it is for MSSM, and it is independent of tan j3. Furthermore, it turns out that the new 
contribution is also insensitive to e'; thus it is nearly the same for cases I and II of the 
fermion mass-matrix. Comparing Eq.(pT|) with Eq.(PBD we see that the new and the standard 
operators are typically quite comparable to one another. Since there is no reason to expect 
near cancellation between them (especially for both v^K^ and Ut-K^ modes), we expect the 
net amplitude (standard-|-new) to be in the range exhibited by either one. It is thus useful 
to obtain the inverse decay rate assuming as if the new operator dominates. Substituting 
Eq.(^T]) into Eq.(^ID and allowing for the presence of the Vt-K^ mode with an estimated 
branching ratio of nearly 0.4 (see Ref. Jl^); one obtains 

V-^{vK+)neu, ~ (1 X lO^Vrs) [8-1/64] {64-l/64}(ir-2 ^ 9 to 1/4). (52) 



The square bracket represents the uncertainty reflected in Eq. (pTD , while the curly bracket 
corresponds to that in the SUSY spectrum and matrix element (Eq. (^OD ). Allowing for a 
net uncertainty at the upper end by as much as a factor of 100 to 600 (say), arising jointly 
from the three brackets in Eq.(^), which can be realized by keeping the SUSY-spectrum 
and the matrix element in the "nearly-central" or "intermediate" range (see below), the new 
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operators related to neutrino masses, by themselves, lead to a proton decay lifetime given 
by: 

r-^(l7ir+)*^^'^^"" ^ (0.7-5) X lO^^rs. (SO(IO) or string G(224))( Indep. of tan/?). (53) 

For instance, taking the curly bracket in Eq.(^) to be 4 to 10 (say) (corresponding for 
example, to /3h = 0.012 GeV^, {my^,/mq) ^ 1/10 to 1/12 and m^- ^ (1 to 1.3)(1.2 TeV)), 
instead of its extreme value of 64, and setting the square bracket in Eq.(0) to be ^ 6, 
and K^"^ ^ 9, which are quite plausible, we obtain: T^^{h'K^)new ~ (2.2 — 4) x 10^^ yrs; 
independently of tan/5, for both MSSM and ESSM. Proton lifetime for other choices of 
parameters, which lead to similar conclusion, are listed in Table 2. 

It should be stressed that the standard d = 5 operators (mediated by the color-triplets 
in the 10^ of SO(IO)) may naturally be absent for a string-derived G(224)-model (see e.g. 
Ref. and 0]), but the new d = 5 operators, related to the Majorana masses of the RH 
neutrinos and the CKM mixings, should very likely be present for such a model, as much as 
for SO(IO). These would induce proton decay 1^. Thus our expectations for the proton decay 
lifetime (as shown in Eq. ([3^jj and the prominence of the n'^K^ mode (see below) hold for 
a string -derived G(224)-model, just as they do for SO(IO). For a string - G(224) - model, 
however, the new d=5 operators would be essentially the sole source of proton decay. 

Nearly the same situation emerges for the case of ESSM embedded in G(224) or SO(IO), 
with low tan/?(Ri 3 to 10, say), especially if the parameters (including (3h, f^w/n^^g, t^qi 
the phase-dependent factor as well as M^fj) happen to be somewhat closer to their extreme 
values so as to extend proton lifetime. In this noted in the previous sub-section, 

the contribution of the standard d=5 operators would be suppressed; and proton decay 
would proceed primarily via the new operators with a lifetime quite naturally in the range 
of 10^^ — 10^^ years, as exhibited above. 



F. The Charged Lepton Decay Modes (p ^'^K^ and p e'^ir^) 

I now note a distinguishing feature of the SO(IO) or the G(224) model presented here. 
Allowing for uncertainties in the way the standard and the new operators can combine with 
each other for the three leading modes i.e. Vt-K~^ , V^K"^ and ^'^K^, we obtain (see Ref. 
for details): 

Bifi+K^) std+new ~ [1% to 50%] K (SO(IO) or string G(224)) (54) 

where k denotes the ratio of the squares of relevant matrix elements for the n'^K^ and VK~^ 
modes. In the absence of a reliable lattice calculation for the i/K~^ mode, one should remain 
open to the possibility of k ~ 1/2 to 1 (say). We find that for a large range of parameters, 



In addition, quantum gravity induced d=5 operators are also expected to be present at some 
level, depending upon the degree of suppression of these operators due to flavor symmetries (see 
e.g. Ref. pi). 
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the branching ratio B{^'^K^) can he in the range of 20 to 40% (if n 1). This prominence of 
the fi~^K^ mode for the SO(10)/G(224) model is primarily due to contributions from the new 
d=5 operators. This contrasts sharply with the minimal SU(5) model, in which the ^^K^ 
mode is expected to have a branching ratio of only about 10~^. In short, prominence of the 
^+^0 j^Q(jg^ if seen, would clearly show the relevance of the new operators, and thereby 



reveal the proposed link between neutrino masses and proton decay [|T9 . 

While the d=5 operators as described here (standard and new) would lead to highly 
suppressed e+vr'' mode, for MSSM or ESSM embedded in SO(IO), the gauge-mediated d=6 
operators, can still give proton decay into e^vr" with an inverse rate ~ XO^^'^^^'^ years, which 
can be as short as about 10^^ yrs. Thus, even within supersymmetric unification, the e+vr'^ 
mode may well be a prominent one, competing favorably with (even) the vK^ mode. 



G. Section Summary 

In summary, our study of proton decay has been carried out within the supersymmetric 
SO(IO) or the G(224)-framework^ with special attention paid to its dependence on fermion 
masses and threshold effects. A representative set of results corresponding to different choices 
of parameters is presented in Tables 1 and 2. The study strongly suggests that, for either 
MSSM or ESSM embedded in SO(IO) or G(224), an upperlimit on proton lifetime is given 
by 

Voton < (1/2-1) X 10^4 yrs, (55) 

with VK^ being the dominant decay mode, and fi~^K^ being prominent. Although there are 
uncertainties in the matrix element, in the SUSY-spectrum, in the phase-dependent factor, 
tan /3 and in certain sensitive elements of the fermion mass matrix, notably e' (see Eq. (^) for 
predictions in cases I versus II), this upper limit is obtained, for the case of MSSM embedded 
in SO(IO), by allowing for a generous range in these parameters and stretching all of them 
in the same direction so as to extend proton lifetime. In this sense, while the predicted 
lifetime spans a wide range, the upper limit quoted above, in fact more like 3 x 10^^ yrs, 
is most conservative, for the case of MSSM (see Eq. (^8]) and Table 1). It is thus tightly 
constrained already by the empirical lower limit on r~^{vK'^) of 1.6 x 10^'^ yrs. For the case 
of ESSM embedded in SO (10), the standard d=5 operators are suppressed compared to the 
case of MSSM; as a result, by themselves they can naturally lead to lifetimes in the range 
of (3 — 10) X 10^^ yrs., for nearly central values of the parameters pertaining to the SUSY- 
spectrum and the matrix element (see Eq.(^)) and Table 1. Including the contribution of 
the new d=5 operators, and allowing for a wide variation of the parameters mentioned above, 
one finds that the range of (10^^ — 10'^^) yrs for proton lifetime is not only very plausible 
but it also provides a rather conservative upper limit, for the case of ESSM embedded in 
either SO(IO) or 0(224) (see Sec. [VIE| and Table 2). Thus our study provides a clear reason 



to expect that the discovery of proton decay should be imminent for the case of ESSM, and 
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As described in Sees. Ill, IV and 
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even more so for that of MSSM. The imphcation of this prediction for a next-generation 
detector is emphasized in the next section. 



VII. CONCLUDING REMARKS 



The preceding sections show that, but for one missing piece - proton decay - the evidence 
in support of grand unification is now strong. It includes: (i) the observed family- structure, 
(ii) the meeting of the gauge couplings, (iii) neutrino-oscillations, (iv) the intricate pattern of 
the masses and mixings of all fermions, including the neutrinos, and (v) the need for B — L 
as a generator, to implement baryogenesis. Taken together, these not only favor grand 
unification but in fact select out a particular route to such unification, based on the ideas 
of supersymmetry, SU(4)-color and left-right symmetry. Thus they point to the relevance of 
an effective string-unified G(224) or SO(10)-symmetry. 

Based on a systematic study of proton decay within the supersymmetric SO(10)/G(224)- 
framework which is clearly favored by the data, and an update as presented here, I have 
argued that a conservative upper limit on the proton lifetime is about (1/2 - l)xlO^^ yrs. 
for the case of either MSSM or ESSM, embedded in SO(IO) or a string - G(224). 

So, unless the fitting of all the pieces listed above is a mere coincidence, and I believe that 
that is highly unlikely, discovery of proton decay should be around the corner. In particular, 
as mentioned in the Introduction, we expect that candidate events should very likely be 
observed in the near future already at SuperK. However, allowing for the possibility that 
proton lifetime may well be near the upper limit or value stated above, a next-generation 
detector providing a net gain in sensitivity by a factor five to ten, compared to SuperK, would 
be needed to produce real events and distinguish them unambiguously from the background. 
Such an improved detector would of course be essential to study the branching ratios of 
certain crucial though (possibly) sub-dominant decay modes such as the fi~^K^ and e~^7i^ as 
mentioned in Sec. |VI F| . 

The reason for pleading for such improved searches is that proton decay would provide 
us with a wealth of knowledge about physics at truly short distances (< 10~^^ cm), which 
cannot be gained by any other means. Specifically, the observation of proton decay, at a 
rate suggested above, with VK~^ mode being dominant, would not only reveal the underlying 
unity of quarks and leptons but also the relevance of supersymmetry. It would also confirm 
a unification of the fundamental forces at a scale of order 2 x 10^^ GeV. Furthermore, 
prominence of the fi'^K'^ mode, if seen, would have even deeper significance, in that in 
addition to supporting the three features mentioned above, it would also reveal the link 
between neutrino masses and proton decay, as discussed in Sec.[V|. In this sense, the role 
of proton decay in probing into physics at the most fundamental level is unique . In view 
of how valuable such a probe would be and the fact that the predicted upper limit on the 
proton lifetime is at most a factor of three to six higher than the empirical lower limit, the 
argument in favor of building an improved detector seems compelling. 

To conclude, the discovery of proton decay would undoubtedly constitute a landmark 
in the history of physics. It would provide the last, missing piece of gauge unification and 
would shed light on how such a unification may be extended to include gravity in the context 
of a deeper theory. 
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APPENDIX: A NATURAL DOUBLET-TRIPLET SPLITTING MECHANISM IN 

SO(IO) 

In supersymmetric SO (10), a natural doublet-triplet splitting can be achieved by coupling 
the adjoint Higgs 45h to a IOh and a 10^, with 45h acquiring a unification-scale VEV in 
the B-L direction [^]: (45h) = (a, a, a, 0, 0) x T2 with a ~ Mu- As discussed in Section |V|, 
to generate CKM mixing for fermions we require (IGh)^ to acquire a VEV of the electroweak 
scale. To ensure accurate gauge coupling unification, the effective low energy theory should 
not contain split multiplets beyond those of MSSM. Thus the MSSM Higgs doublets must be 
linear combinations of the SU(2)£, doublets in IOh and 16h- A simple set of superpotential 
terms that ensures this and incorporates doublet-triplet splitting is [|I^]: 

Wh = \ 10h45h 10'h + Mio 10'^^ + X'TQhI^h IOh + Mie I^hI^h • (Al) 

A complete superpotential for 45h, 16h, 16h, 10^, 10'^;^ and possibly other fields, which 
ensure that (a) 45h, 16h and 16h acquire unification scale VEVs with (45h) being along the 
{B-L) direction; (b) that exactly two Higgs doublets {Hu, H^) remain light, with being a 
linear combination of (IOh)^ and (IGh)^; and (c) there are no unwanted pseudoGoldstone 
bosons, can be constructed. With (45h) in the B-L direction, it does not contribute to 
the Higgs doublet mass matrix, so one pair of Higgs doublet remains light, while all triplets 
acquire unification scale masses. The light MSSM Higgs doublets are 



= lOu, Hd = cos7l0rf + sin7l6rf, (A2) 

withtan7 = A'(16h)/Mi6. Consequently, (10)^ = (cos7)frf, (16^) = (sin7)t>^, with {Hd) = 
Vd and (16^) and (10^) denoting the electroweak VEVs of those multiplets. Note that 
is purely in IOh and that (10^)^ + (16^)^ = fj. This mechanism of doublet-triplet (DT) 
splitting is the simplest for the minimal Higgs systems. It has the advantage that meets 
the requirements of both D-T splitting and CKM-mixing. In turn, it has three special 
consequences: 

(i) It modifies the familiar SO(10)-relation tan/3 = Vu/vd = rrit/mh ~ 60 to 0: 



It is worth noting that the simple relationship between cos 7 and tan/3 - i.e. cos 7 ~ 
tan P / {rrit / mfj) - would be modified if the superpotential contains an additional term like \"16h ■ 
16h ■ lOj:/, which would induce a mixing between the doublets in 10^, 16^ and 10^. That in turn will 
mean that the upper limit on M(,ff cos 7 following from considerations of threshold corrections (see 
below) will not be strictly proportional to tan /?. I thank Kaladi Babu for making this observation. 
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tan/?/ COS7 ~ w^t/'o^fe ~ 60 (A3) 

As a result, even low to moderate values of tan/5 ^ 3 to 10 (say) are perfectly allowed in 
SO(IO) (corresponding to C0S7 ~ 1/20 to 1/6). 

(ii) The most important consequence of the DT-splitting mechanism outlined above is 
this: In contrast to SU(5), for which the strengths of the standard d=5 operators are pro- 
portional to (MhJ"^ (where ~ few x 10^^ GeV (see Eq. (H)), for the SO(10)-model, 
they become proportional to M~j^j, where M^jf = (Aa)^/Mio' ~ M'^/Miq^. As noted in 
Ref. [Q, MiQ! can be naturally smaller (due to flavor symmetries) than Mx and thus M^ff 



correspondingly larger than Mx by even one to three orders of magnitude. Now the proton 
decay amplitudes for SO (10) in fact possess an intrinsic enhancement compared to those for 
SU(5), owing primarily due to differences in their Yukawa couplings for the up sector (see 
Appendix C in Ref. [T^). As a result, these larger values of Mg// ~ (10^^ — 10^^) GeV are 
in fact needed for the SO(10)-model to be compatible with the observed limit on the proton 
lifetime. At the same time, being bounded above by considerations of threshold effects (see 
below), they allow optimism as regards future observation of proton decay. 

(iii) Mgff gets bounded above by considerations of coupling unification and GUT-scale 
threshold effects as follows. Let us recall that in the absence of unification-scale threshold 
and Planck-scale effects, the MSSM value of a^imz) in the MS scheme, obtained by assuming 
gauge coupling unification, is given by a3(^z)MSSM — 

0.125 - 0.13 [0]. This is about 5 to 
8% higher than the observed value: as^mz) = 0.118 ± 0.003 [jl3|. Now, assuming coupling 
unification, the net (observed) value of 0:3, for the case of MSSM embedded in SU(5) or 
SO(IO), is given by: 

asimzUt = «3(mz)MSSM + Aasimz)^^^"^ + A'3 (A4) 

where Aa3(m^)DT and A3 represent GUT-scale threshold corrections respectively due 
to doublet-triplet splitting and the splittings in the other multiplets (like the gauge 
and the Higgs multiplets), all of which are evaluated at m^. Now, owing to mix- 
ing between 10^ and 16^ (see Eq. (|A2|) ), one finds that Aa3{mz)DT is given by 
(a3(mz) V2vr) (9/7) ln(Meff cos j/Mx) H . 



As mentioned above, constraint from proton lifetime sets a lower limit on M^s given by 
Meff > (1 - 6) X lO^^GeV. Thus, even for small tan/? 2 (i.e. cos 7 ^ tan(/3/60) ^ 1/30), 
Aa;3(mz)DT is positive; and it increases logarithmically with Mes- Since Oi^(fnz)'^<^^y[ is 
higher than a3{mz)ohs, and as we saw, Aa3{mz)DT is positive, it follows that the corrections 
due to other multiplets denoted by 6^ = Aya^^mz) should be appropriately negative so 
that a2,{mz)nct would agree with the observed value. 

In order that coupling unification may be regarded as a natural prediction of SUSY uni- 
fication, as opposed to being a mere coincidence, it is important that the magnitude of the 
net other threshold corrections, denoted by ^3, be negative but not any more than about 8 



to 10% in magnitude (i.e. — ^3 ^ (8 — 10)%). It was shown in Ref. fl^ that the contributions 
from the gauge and the minimal set of Higgs multiplets (i.e. 45/^, 16^, 16h and 10^) leads 
to threshold correction, denoted by 6'^, which has in fact a negative sign and quite naturally 
a magnitude of 4 to 8%, as needed to account for the observed coupling unification. The cor- 
rection to a^{mz) due to Planck scale physics through the effective operator Fh^F^'^A5h/M 
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does not alter the estimate of because it vanishes due to antisymmetry in the SO(IO)- 
contraction. 

Imposing that ^3 (evaluated at mz)he negative and not any more than about 10-11% 
in magnitude in turn provides a restriction on how big the correction due to doublet-triplet 
splitting - i.e. Aa3{mz)r,T - can be. That in turn sets an upper limit on M^jf cos'j, and 
thereby on Me// for a given tan (3. For instance, for MSSM, with tan /3 = (2, 3, 8), one obtains 
(see Ref. fll): M^// < (4,2.66, 1) x lO^^GeV. Thus, conservatively, taking tan/? > 3, one 
obtains: 

Meff < 2.7 X lO^^GeV (MSSM). (A5) 



Limit on Me// For The case of ESSM 

Next consider the restriction on Me// that would arise for the case of the extended su- 
persymmetric standard model (ESSM), which introduces an extra pair of vector-like families 
(16 + 16) of SO(IO)) at the TeV scale pO[(see also footnote 11). In this case, a^nif is raised 
to 0.25 to 0.3, compared to 0.04 in MSSM. Owing to increased two-loop effects the scale 
of unification Mx is raised to (1/2 — 2) x lO^^GeV, while a3('^z)EssM lowered to about 



0.112-0.118 |20|J5^. 



With raised Mx, the product Me//cos7 ~ Me//(tan/3)/60 can be higher by almost a 
factor of five compared to that for MSSM, without altering Aa3{mz)-DT- Furthermore, since 
Q^3("^z)mssm is typically lower than the observed value of asi^rnz) (contrast this with the case 
of ESSM), for ESSM, Me// can be higher than that for MSSM by as much as a factor of 2 
to 3, without requiring an enhancement of ^3. The net result is that for ESSM embedded in 
SO(IO), tan/3 can span a wide range from 3 to even 30 (say) and simultaneously the upper 
limit on Me// can vary over the range (60 to 6)xlO^^Ge\^, satisfying 

Me// < (6 X 10^^GeV)(30/tan/5) (ESSM), (A6) 

with the unification-scale threshold corrections from "other" sources denoted by ^3 = 
A3/a3(m^) being negative, but no more than about 5% in magnitude. As noted above, 
such values of ^3 emerge quite naturally for the minimal Higgs system. Thus, one important 
consequence of ESSM is that by allowing for larger values of Me// (compared to MSSM), 
without entailing larger values of ^3, it can be perfectly compatible with the limit on proton 
lifetime for almost central values of the parameters pertaining to the SUSY spectrum and 
the relevant matrix elements (see Eq.(40)). Further, larger values of tan/3 (10 to 30, say) 
can be compatible with proton lifetime only for the case of ESSM, but not for MSSM. These 
features are discussed in the text, and also exhibited in Table 2. 
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TABLE 2. VALUES OF PROTON LIFETIME {T-\p 

RANGE OF PARAMETERS 



uK+)) FOR A WIDE 



Parameters 


MSSM - 
Std. 


SO(IO) 
d— 5 


ESSM ^ SO(IO) 
Std. d=5 


J MSSM' 
1 ESSM 


> ^ G(224)/SO(10) 


eiGiiiGiii J 


iiitei meQ. 


t 06 pnase 


Intermed. e' & phased 


New d=5tt 




tan/3=3 


tan/3=20 


tan /3=5 


tan /5=20 


Independent of tan (3 


Nearly "central" 


0.7 X 10^2 


1.6x1030 


1.1x1034 


0.7x1033 


0.7 X 1033 


{ }=2 


yrs 


yrs 


yrs* 


yrs 




yrs^T 


Intermediate 


2.8x1032 


0.6x1031 


0.4x1035 


2.8x1033 




2.8x1033 


{ }=8 


yrs 


yrs 


yrs* 


yrs 




yrs''"'' 


Nearly Extreme 


1.1X1033 


2.6x1031 


1.7x1035 


1.1X1034 




1.1x1034 


{ }=32 


yrs 


yrs 


yrs* 


yrs 




yrs^''" 



"In this case, lifetime is given by the last column. 

• Since we are interested in exhibiting expected proton lifetime near the upper end, we 
are not showing entries corresponding to values of the parameters for the SUSY spectrum 
and the matrix element (see Eq. f^OD , for which the curly bracket appearing in Eqs. f^Tf) , (|49|) , 
(0)) would be less than one (see however Table 1). In this context, we have chosen here 
"nearly central", "intermediate" and "nearly extreme" values of the parameters such that 
the said curly bracket is given by 2, 8 and 32 respectively, instead of its extreme upper-end 
value of 64. For instance, the curly bracket would be 2 if (3h = (0.0117) GeV3, ^ 1.2 

1.44 TeV and 
i ^2(1.2 TeV) 



ma 



TeV and m^/nig ^ (1/7.2), while it would be 8 if Ph = 0.010 GeV3, 
m^/nig ^ 1/10; and it would be 32 if, for example, Ph = 0.007 GeV3, 
and m^/rriq ~ 1/12. 

I All the entries for the standard d=5 operators correspond to taking an intermediate 
value of e' ~ (1 to 1.4) x 10^^ opposed to the extreme values of 2 x 10"^ and zero for cases 
I and II, see Eg. (|3^) ) and an intermediate phase-dependent factor such that the uncertainty 
factor in the square bracket appearing in Eqs.(^) and (49) is given by 5, instead of its 
extreme values of 2 x 4 = 8 and 2.5 x 4 = 10, respectively. 

If For the new operators, the factor [8-1/64] appearing in Eq.(^) is taken to be 6, and 
K~^, defined in SecjV^A, is taken to be 9, which are quite plausible, in so far as we wish to 
obtain reasonable values for proton lifetime at the upper end. 

• The standard d=5 operators for both MSSM and ESSM are evaluated by taking the 
upper limit on Me// (defined in the text) that is allowed by the requirement of natural 
coupling unification. This requirement restricts threshold corrections and thereby sets an 
upper limit on Mg//, for a given tan/3 (see Sec.[V3 and Appendix). 

• For all cases, the standard and the new d=5 operators must be combined to obtain the 
net amplitude. For the three cases of ESSM marked with an asterisk, and other similar cases 
which arise for low tan/3 fa 3 to 6 (say), the standard d=5 operators by themselves would 
lead to proton lifetimes typically exceeding (0.1 — 0.7) x 103^ years. For these cases, however, 
the contribution from the new d=5 operators would dominate, which quite naturally lead to 
lifetimes in the range of (1033 — 103^) years (see last column). 

• As shown above, the case of MSSM embedded in SO (10) is tightly constrained by 
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present empirical lower limit on proton lifetime (Eg. P^)). In this case, only low values of 
tan/3 < 3, with the parameters (pertaining to the SUSY spectrum, matrix element and 
phase-dependent factor) having their "nearly extreme" or extreme values (as in Eq.(POD) 
can lead to lifetimes in the range of (1 — 3) x lO^^yrs (see Table and Eg . (|48|) ) , compatible 
with present empirical limit. Other cases of MSSM - especially with tan/3 > 5 and/or 
"nearly central" or even "intermediate" range of parameters - seem to be excluded, subject 
(of course) to our reguirement for natural coupling unification (see Sec.0 and Appendix). 

• Including contributions from the standard and the new operators, the case of ESSM, 
embedded in either G(224) or SO(IO), is, however, fully consistent with present limits on 
proton lifetime for a wide range of parameters; at the same time it provides optimism that 
proton decay will be discovered in the near future, with a lifetime < 10^^ years. 

• The lower limits on proton lifetime are not exhibited. In the presence of the new opera- 
tors, these can typically be as low as about 10^^ years (even for the case of ESSM embedded 
in SO(IO)). Such limits and even higher are of course long excluded by experiments. 

• Allowing for a wide variation in the relevant parameters, we thus see that a conservative 
upper limit on proton lifetime is given by the range of (1/2 — 1) x 10^^ years for ESSM and 
(of course) MSSM, embedded in SO (10) or string-G(224). 
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